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INTTRODUCriON 


The  main  theme  of  the  present  Grant  was  the  noninvasive  spectral  analysis  of 
ultrasonic  fields  where  nonlinear  properties  of  the  medium  of  propagation  cause 
initially  sinusoidal  ultrasonic  waves  to  become  distorted  and  thus  exhibit  more 
than  one  frequency  component  in  the  sound  field. 

There  are  a  number  of  cases  where  such  "finite  amplitude  effects"  -  as 
opposed  to  infinitesimally  small  amplitudes  of  the  signal  where  nonlinearities 
can  be  neglected-  occur.  All  of  them  are  related  to  high  amplitude  waves  propa¬ 
gating  in  water;  the  most  common  ones  are: 

High  amplitude  continuous  wave  bounded  beams, 

Ultrasonic  pulses. 

Focused  ultrasonic  beams. 

The  work  performed  consisted  primarily  of  theoretical  and  experimental 
analyses  of  the  spectral  contents  of  these  three  sound  fields  and,  secondarily, 
of  the  investigation  of  reflection,  transmission,  and  reradiation  phenomena  when 
such  fields  encounter  liquid/solid  interfaces,  like  solid  plates  immersed  in 
water. 


The  theoretical  analysis  was  concerned  with  two  topics.  The  first  was  the 
computation  model  of  the  changes  in  the  harmonic  contents  and  the  pressure 
profile  of  the  sound  signal  as  a  function  of  travel  distance. 

(a)  In  the  case  of  line  focused  ultrasotmd  the  theoretical  sound  pressure 
field  predictions  were  based  on  a  numerical  solution  of  a  differential  equation 
describing  nonlinear  effects,  the  Khokhlov-Zabolotskaya-Kuznetsov  equation.  A 
method  was  found  to  greatly  enhance  the  cc’::putational  effectiveness  of  the 
numerical  solution. 

(b)  In  the  case  of  pulsed  ultrasound,  the  important  factors  were  f(»ind  to 
be  the  Fourier  series,  whose  lowest  frequency  component  is  the  pulse  repetition 
rate,  and  the  changes  in  that  frequency  spectrum  due  to  nonlinear  processes  as 
the  signal  progresses. 


The  second  general  theoretical  topic  was  concerned  with  the  calculation  of 
the  acousto-optic  interact  ion  of  pulsed  ultrasound  and  monochrooatic  light.  This 
interaction  produces  a  distinctive  diffraction  pattern  where  the  light  intensity 
in  the  many  diffraction  orders  is  determined  by  the  frequency  spectrum  of  the 
signal.  The  various  aspects  of  these  two  theoretical  approaches  will  be  dis¬ 
cussed  in  later  sections  of  this  report. 

The  experimental  analysis  was  done  by  using  acousto-optic  techniques  to 
probe  the  sound  fields.  The  advantage  here  is  the  fact  that  the  sound  field  is 
not  being  disturbed  by  the  measurement  since  no  material  object,  like  a  hydro¬ 
phone,  needs  to  be  placed  in  the  sound  field.  Roth  diffraction  order  analysis 
and  Schlieren  techniques  were  used  to  verify  theoretical  predictions  and  in  some 
cases  to  establish  the  basis  for  further  theoretical  considerations. 

This  report  is  structured  according  to  the  outline  given  above.  The  sec¬ 
tions  are  separated  and  are  preceded  by  individual  Summary  Pages .  The  Appendices 
contain  listings  of  computational  codes  and/or  reprints  of  papers  which  have  a 
bearing  on  the  main  themes. 
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This  section  was  prepared  by  Dr.  Thomas  Neighbors  who  has  done  the  computational 
approach  to  the  calculation  of  pressure  profile  in  a  line-focused  ultrasonic 
beam.  The  experimental  results  presented  in  this  section  were  obtained  by 
Herbert  Ruf.  Both  these  authors  were  partially  supported  by  the  present  Grant; 
Dr.  Neighbors  for  a  very  short  time  as  Research  Associate  Professor,  and  Dr.  Ruf 
for  most  of  the  duration  of  the  Grant  as  graduate  student  and  later  as  post¬ 
doctoral  fellow. 

This  section  is  entitled  ’’Line  Focused  Transducers  -  Equations  and  Field 
Determinations.”  The  KZK  equation  is  stated  and  then  used  to  apply  to  line- 
focused  sound  sources.  The  needed  mathematical  operations  are  indicated  so  the 
equation  can  be  written  in  normalized  terms. 

After  the  focused  beam  boundary  conditions  are  enumerated,  the  KZK  equation 
Fourier  decomposition  for  line  focused  sources  is  performed.  The  main  problem 
is  the  solution  of  an  infinite  set  of  coupled  nonlinear  equations.  The  key  steps 
to  do  this  in  an  efficient  yet  reliable  way  are  listed  and  the  numerical  solu¬ 
tions  are  discussed,  concentrating  on  the  validity  of  the  idea  of  accepting 
"floating"  boundary  conditions,  the  flexible  allocation  of  harmonic  terms,  and 
lastly,  the  convergence  of  the  solution. 

This  is  followed  by  the  listing  of  results  for  concrete  examples  of  line- 
focused  beams,  giving  the  initial  conditions  at  the  source.  This  is  compared 
with  experimental  results,  where  the  experiments  are  performed  making  use  of 
acousto-optic  methods. 

This  section  closes  with  a  statement  of  the  predicted  changes  in  the  focal 
length  of  the  transducer  when  the  pressure  is  increased. 
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Line  Focus  Equation  Derivation 


1.0  Introductk>n 

The  Khokhlov-Zabolotskaya-Kuznetsov  (KZK)  parabolic  equation  is  generally  used 
to  account  for  nonlinearity,  diffraction  and  attenuation  in  a  directive,  finite-amplitude  sound 
beam*.  The  KZK  equation  is  given  by: 


b  ci^P_  P 

*ar'  2  ‘ 


where 

Cq  =  isentropic  speed  of  sound, 

Po  =  fluid  density, 

P  =  coefficient  of  nonlinearity, 

t  =  retarded  time  (t  -  z/c^), 

t  =  time, 

z  =  on-axis  distance, 

P  =  pressure, 

=  transverse  Laplican,  and 
b  =  thermo-viscous  loss  term. 

Originally  Eq.  1  was  derived  without  thermo-viscous  losses  by  Khokhlov  and 
Zabolotskaya.  The  adaption  to  include  dissipation  is  due  to  Kuznetsov.  A  recent  re- 
derivation^  has  shown  that  Eq.  1  accounts  for  nonlinearity,  diffraction  and  absorption  to  an 
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equal  order  of  magnitude.  The  term  corresponds  to  dif&action  effects.  In  deriving 
Eq.  1  a  well  collimated  beam  is  assumed,  e.g.  ka  >>  1  with  k  =  cu/Cg  and  with  a  ==  the 
source  radius. 

The  following  sections  outline  the  approach  used  in  applying  the  KZK  equation  to 
a  line  focused  radiator.  Sect'on  2  starts  by  outlining  the  modifications  in  the  parabolic 
equation  necessary  for  its  initial  application  to  the  line  focused  radiator.  Section  3  briefly 
outlines  the  change  in  the  boundary  conditions  necessary  when  the  KZK  equation  is  used 
for  a  focused  source.  This  is  followed  in  section  4  by  the  Fourier  series  expansion  of  the 
modified  KZK  equation  which  results  in  an  infinite  set  of  coupled  nonlinear  differential 
equations  for  the  pressure  spectral  terms.  Section  5  takes  the  coupled  nonhnear  differential 
equations  and  discusses  the  computational  efficiency  enhancement  techniques  which  allow 
their  numerical  solution  on  a  desk  top  computer. 

2.0  KZK  Equation  For  Line  Focused  Source 

For  a  geometry  which  considers  only  the  y  z  plane  with  z  the  propagation  direction 
and  y  perpendicular  to  z,  the  transverse  gradient  can  be  replaced  with  d‘y.  This  approach 
assumes  that  to  first  order  the  variations  in  the  x  direction  can  be  ignored. 


dzdt'  2dy^  2poCo  at'* 


If  the  retarded  time  is  replaced  by  a  normalized  retarded  time,  i.e.  r  =  mt*,  then  dt* 
=  drlcD  and  equation  (2)  can  be  rewritten  as 


2 


(3) 


t>u>*  e?p_ 

dzdx  2  dy^  2p^o  2poCo 


The  geometry  used  in  applying  the  KZK  equation  to  focused  sound  propagation  is 
in  Figure  1.  The  acoustic  field  originates  at  a  concave  source,  of  transverse  diameter  2a, 
and  propagates  in  the  z  direction  towards  the  geometric  focus  w^kh  is  nominally  located 
at  the  source  radius  of  curvature  d. 


To  simplify  equation  (3),  renormalize  the  on-axis  distance,  z,  to  the  radius  of 
curvature  of  the  line  focus  radiator,  d,  with  the  new  distance  parameter  being  a  =  z/d. 


<a  ^P  Cq^P  ^P  _  ^P^ 

d  dodx  2  dy^  2poCo  ^P(fo 

Equation  (4)  can  be  reduced  even  further  by  normalcing  the  transverse  coordinate 
to  the  radiator’s  transverse  radius  a.  as  shown  in  Rgure  1.  With  the  new  transverse 
coordinate  r  =  y/a,  equation  (4)  becomes. 
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(5) 


CJ  gp  ^P_ 

'f  dodx  2a*  A-*  2poCo 


At  this  point,  we  can  convert  to  a  normalized  pressure  by  dividing  P  by  the  on-axis 
peak  pressure  at  the  surface  of  the  line  focused  radiator,  f\„  e.g.  p  =  P/p^ 


“Pq  c<Pq^p 

d  dadx  2a*  a-*  2p^cldt^  Ip^cl  dr* 


Multiply  by  d/rop,,  and  we  have 


df^p_feto*d^p.P<*>Po^d*p* 
dear  2«a*ar*  2poCo^*  2poCo 


Equation  (7)  can  be  simplified  by  using  the  standard  definitions  of  the  Rayleigh 
distance,  Zq,  and  the  discontinuity  distanced^  given  in  equations  (8)  and  (9). 


and 


j  .  pQgQ 

'  PwRo 


(«) 


(9) 


Equation  (10)  is  the  result  of  the  substitution. 
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(10) 


d  _  d 

dadv  4zq  dr^  2pf^ci  3t*  2/^  dr^ 

Identifying  the  last  set  of  constants  as  a,  the  thermal  viscous  loss  coefficient  and 
making  the  substitution  shown  below  for  ad  yields  equation  (12). 


ad= 


2p<jCo 


Then  we  have  the  line  focus  KZK  equation  in  normalized  terms. 


(11) 


^P _ d  ^P  ^ 

do dr  42q  dr^  dz^  2/^ 


(12) 


3.0  Focused  Beam  Boundary  Conditions 

To  effect  a  solution  of  Eq.  12  it  is  necessary  to  specify  the  initial  pressure  distribution 
p(y,z,t)  at  z  =  0.  The  standard  approach  for  specifying  the  initial  condition  uses  the  quasi- 
optical  approximaticHi  which  modifies  the  amplitude  and  phase  of  the  pressure  distributimi 
in  the  xy  plane  at  z  =  0  to  produce  the  amplitude  and  phase  distribution  that  would  occur 
at  the  surface  of  the  concave  source.  Figures  2a  and  2b  flhistrate  the  phase  and  amplitude 
adjustment,  respectively.  The  phase  adjustment  shown  in  Figure  2a  is  k*dx  where  k  =  InIX 
and  dx  is  the  distance  that  would  have  to  be  propagated  from  the  xy  plane  to  the  curved 
surface. 

The  amplitude  correction  requires  a  bit  more  explanation  as  shown  in  Fig  2b.  Since 
the  initial  pressure  distribution  in  each  elementary  radial  segment  along  the  y  axis,  dy,  over 
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the  interva]  (0,a]  represents  the  pressure  on  a  segment  of  the  curved  surface,  the  required 
increase  in  the  normalized  pressure  at  z  =  0  is  given  by  D  •  dd/dy.  As  shown  in  Figs.  2a 
and  2b,  0  is  the  angle  measured  from  the  z  axis  and  D  is  the  source  radius  of  curvature. 


(•)  Phase  Correction  (b)  AmpUinde  Correction 

Figure  2.  Initial  Pressure  Phase  and  Amplitude  Adjustment 

4.0  KZK  equation  Fourier  decomposition  for  line  focused  source 

Since  a  general  analytic  solution  of  the  KZK  equation  has  not  been  found  to  date^ 
the  equation  is  usually  solved  by  a  number  of  approximatimis  and  numerical  techniques. 
The  approach  used  here  has  been  adapted  from  the  numerical  sohiticm  initially  developed 
by  Aanonsen*,  enhanced  Hart  and  Hamilton’,  and  revised  by  Ne^hbots  and  Bjomo. 
The  numerical  approach  allows  solutions  within  a  range  of  assumptions  which  are  not  too 
restrictive. 

Assume  the  initial  pressure  at  <7  =  0  can  be  represented  as  a  periodic  functicm  of  the 


(13) 


with 


fp(^)e-^dx 


(14) 


Now  if  we  substitute  Eq.  13  into  Eq.  12  we  have  the  following  set  of  coupled 
nonlinear  differential  equations  in  terms  of  the  Fourier  components: 


.4>.  d  a’p. 


do 
_d 

21,  dx^ 


ad(/n)V>, 


4Zo  dr^ 


(15) 


To  convert  Eq  IS  into  a  set  of  coupled  partial  differential  equations  in  terms  of  p, 
requires  a  transformation  of  the  right  hand  side  of  Eq  15.  First,  let  k  -  n  +  m,  then  we 
have 


E  i  Pm^Pm^^i  T,PJ*k-n^ 


Rotating  the  indices  and  letting  k  -»  n  and  n  ^  k  yields 


The  products  of  e^^  now  become  an  infinite  set  of  coupled  partial  differential 


jlPlPn-k^^ 

*— 


(17) 
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equations  for  p,  of  the  form 


.  d  ^Pm 

jn — = - - 

do  4zq  dr^ 


-ad(jn)\= 


dpnf 

2/. 


EPiP.^ 


(18) 


Since  the  initial  pressure  is  real,  p,  =  p*^  where  the  *  denotes  the  complex 
conjugate.  Using  this  property  and  manipulating  the  summation  indices  yields  equation  (19) 
uliich  is  an  infinite  set  of  coupled  partial  differential  equations  for  the  Fourier  ampbtude 
coefficients  p,. 


.  d  ^Pn 

jn — =- - - 

do  4z^  dr^ 


ad(jnfp^= 


dijnf 


21. 


"•-1  •• 

EPiP-t*^  E  PiPk» 


(19) 


The  standard  approach  for  solving  the  infinite  set  of  coupled  nonlinear  equations 
truncates  the  infinite  spectral  series,  hneaiizes  the  equations,  and  applies  a  numerical  finite 
difference  scheme.  The  key  steps  in  this  solution  are  discussed  next. 


As  a  first  step,  write  p,  as  its  real  and  imaginary  parts, 

Pm=K-Kn 

where 

S 

A,=^//»(T)COS(»T>fT 


(20) 


(21) 


and 
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(22) 


a 

^,=-j‘p(t)sm(nT>/T. 


Substituting  Eq.  22  into  Eq.  19,  and  separating  real  and  imaginary  parts  yields  the 
infinite  set  of  coupled  partial  differentia]  equations  for  h,  and  gg  given  in  Eq.  23  and  24. 


+eLdn^h. 


21 


4/iZo  dr^ 


dn 


Jb>l 


k-**l 


(23) 


da 


d 

4/tZo  dr  2 


+adn\ 


a-l 


k-n^l 


(24) 


Given  an  axi-symmetric  source,  wHb  time  history  such  as  sinr,  the  coupled  differential 
equations  describe  the  growth  of  harmonics  due  to  nonlinearity  and  the  modification  o{  the 
sotmd  field  due  to  diffraction  effects.  Since  the  fundamental  frequency  corresponds  to  n 
=  1,  the  attenuation  of  these  harmonics  is  determined  by  om^  which  is  equivalent  to  a 
square  law  attenuation  term.  The  solutions  can  be  adjusted  to  account  for  ntm-square  law 
attenuation  terms  by  replacing  n^  by  n*  where  g  is  the  attenuaticm  power  law. 

5.0  KZK  Equation  Nmnerical  Solution 

As  discussed  earlier,  since  no  closed  form  solution  exits  for  the  KZK  equation  the 
usual  approach  k  a  numeric^  solution.  The  general  approach  corrsists  of  transforming  Eqs. 
23  and  24  into  a  series  of  linearized  difference  equations  which  can  be  solved  using  standard 
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numerical  techniques^  Since  ibe  specific  numerical  solution  is  discussed  in  detail  in  a 
number  of  references,  this  section  first  focuses  on  the  unique  enhancements  in 
computational  efficiency  that  are  implemented  in  current  work  -  the  dynamic  allocation  Of 
harmonic  terms  in  during  the  computation  and  the  use  of  a  floating  boundary  condition 
which  increases  computational  throughput  whfle  avoiding  some  of  the  artificial  boundary 
problems  discussed  at  the  12'*'  ISNA^  and  second  discusses  the  convergence  of  the 
numerical  solution. 

5.1  "Floating"  or  Expanding  Boundary  Condition 

A  traditional  problem  in  finite  difference  calculations  is  the  initial  selection  of  the 
outer  mesh  size  for  the  calculation.  This  is  important  since  the  pressure  front  eiqrands 
radially  as  the  sound  propagates.  The  improper  selection  of  the  outer  boundary  of  the 
calculational  mesh  can  result  in  one  of  two  of  the  possible  outcomes: 

•  The  creation  of  an  artificial  rigid  boundary  ^ich  creates  an 
spurious  reflections  as  the  signal  propagates.  Ihese  reflections 
in  turn  "poison"  the  answer. 

•  An  excessively  slow  calculation  due  to  a  large  mesh  size.  The 
calculation  rate  is  inverse^  proportional  to  the  number  of 
transverse  mesh  pants  in  the  finite  difference  calculation. 

A  simplistic  approach  can  be  used  to  prevent  the  inclusion  of  an  artificial  boundary 
and  at  improve  computational  efficiency.  This  aj^roach  is  shown  symbolical^  in  Figure  3. 
When  the  pressure  fundamental  frequency  qjectral  amplitude  magnitude  at  the  third  point 
before  the  edge  of  the  boundary,  A[n-3],  exceeds  a  pre-specified  threshold  value  two 
additional  computational  cells  are  added  to  the  calculation.  The  expanding  pressure  front 
does  not  reach  the  rigid  edge  of  the  boundary.  By  expanding  the  boundary  when  required 
the  calculation  time  can  be  reduced  by  as  much  as  a  factor  of  two,  depending  on  the  problem. 
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Figure  3.  "Floating"  boundaiy  condition  numerical  technique 

An  example  of  the  effectiveness  of  this  technique  is  shown  in  Figure  4.  In  this 
instanced,  a  line  focused  transducer  with  a  30°  aperture  angle  is  located  at  z  =  0.  The 
boundary  at  the  start  of  the  calculations  was  inhiaUy  set  to  be  10%  larger  than  the 
transducer  physical  boundary  and  the  calculational  threshold  was  varied  between  l.OxlO"* 
and  1.0x10  '^.  For  the  lowest  threshold  value  the  computational  boundaiy  expands  by  about 
a  factor  of  10  when  propagating  from  the  face  of  the  transducer  to  40%  beyond  the 
geometric  focus.  As  the  threshold  for  expanding  the  boundaiy  is  reduced  the  boundaiy 
expands  at  a  faster  rate,  e.g.  for  a  threshold  of  1.0x10'*^  the  boundaiy  has  expanded  an  order 
of  magnitude  at  a  normalized  distance  70%  beyond  the  face  of  the  transducer.  Allowing  the 
boundaiy  condition  to  "float’  increases  the  computational  efficiency  by  a  factor  two.  A 
relatively  simple  and  effective  technique  for  improving  computational  efficieniy. 

A  natural  question  is  "how  mudi  error  is  introduced  by  this  approximation?”  This 
can  be  addressed  by  examim'ng  the  noimalized  pressure  amplitude  for  the  propagating 
sound  field  over  the  range  from  z  =  0.4  to  1.4  at  a  transverse  location  lined  up  with  the 
edge  of  the  transducer,  e.g  1.0  in  Figure  4.  These  results  are  shown  in  Table  1  for  four 
threshold  values  -  l.OxlO"*,  1.0x10^,  1.0x10^,  and  LOxlO"*®.  Once  the  threshold  is  deceased 
to  ~  1.0x10'^  there  is  essentially  no  difference  in  the  numerical  result 
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Figure  4.  "Floating"  boundaiy  condition  expansion 


Table  1  >  Boundaiy  Threshold  Impact  On  Precision 


z 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 


10“* 

0.036769598 

0.027677476 

0.038064044 

0.043991060 

0.011198936 

0.039476649 

0.063521360 

0.057549026 

0.030258185 

0.073250350 

0.110690431 


THRESHOLD 

1.0-* 


0.038769408 

0.027674593 

0.038083162 

0.043M5829 

0.011194064 

0.039474818 

0.063523095 

0.057547647 

0.030259337 

0.073248796 

0.110664937 


I.O"* 
0.038769408 
0.027674592 
0.038083161 
0.043095831 


I.O** 

0.038769408 

0.027674592 

0.038083161 

0.043995831 


0.01119406S 

4>.039474820 

0.063523093 

0.057547643 

0.030259324 

0.073248786 

0.110684927 


0.011194065 

0.039474820 

0.063523003 

0.057547643 

0.030259324 

0.073248786 

0.110684927 


12 


5.2  Dynamic  AlkKation  Of  Hannonic  Terms 

The  classic  problem  encountered  in  the  numerical  sohiticm  of  ocmlinear  equations 
using  a  Fourier  representation  is  the  large  number  of  harmonic  terms  that  are  required  due 
to  harmonic  growth.  If  two  few  terms  are  used  in  the  Fourier  series  which  represents  the 
solution,  a  truncation  error  occurs  in  last  few  terms  which  start  to  grow  in  magnitude 
abnormally.  As  this  truncation  error  continues  to  increase  it  begins  to  affect  the  preceding 
terms.  Eventually,  the  truncation  error  distorts  the  solution.  Generally,  to  avoid  this 
problem  a  number  of  alternatives  are  available,  such  as: 

•  increase  the  number  of  harmonics  to  push  the  truncation  error 
further  out, 

•  size  the  number  of  harmonics  to  the  precision  required  and 
increase  the  attenuation  of  higher  harmonics  to  compensate,  or 

•  use  either  of  the  preceding  and  filter  the  output  by  discarding 
a  number  of  higher  harmonic  terms  where  the  truncation  error 
is  obvious. 

In  using  any  of  these  approaches,  a  sufiBcient  number  of  haimmics  has  to  be 
retained  to  provide  the  basic  precision  required  in  the  output.  In  general  regardless  of  the 
approach  used,  the  computational  eCGciency  is  proportional  to  the  number  of  harmonics 
squared. 

A  simple  computational  sohiti<xi  is  available  cons^ent  with  the  physics  oi  the 
process.  The  basic  algorithm  is  shown  symbolically  in  Figure  3.  At  the  beginning  of  the 
computation  a  minimum  number  of  harmonics  are  allowed  in  the  computation.  All 
harmonic  terms  are  initialized  to  zero.  At  each  longitudinal  step  in  the  computation,  the 
magnitude  of  the  last  harmonic  is  checked,  if  it  is  above  a  preset  threshold  another 
harmonic  term  is  added. 
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Figure  3.  Dynamic  harmcmic  allocation  process. 

The  threshold  is  selected  as  a  fraction  of  the  normalized  peak  initial  pressure 
amplitude.  The  impact  of  different  threshold  choices  is  shown  in  Tables  Z  3,  and  4.  At  the 
beginning  of  the  calculation  the  maximum  number  of  harmonic  terms  is  set  to  four.  Table 
3  is  the  on-axis  amplitude  spectra  at  z  =  0.4  when  the  threshold  is  varied  from  1.0x10^  to 
1.0x10  '^.  The  amplitude  spectra  magnitude  displayed  to  four  sigm'ficant  figures.  N  is  the 
number  of  the  spectral  term  with  1  corresponding  to  the  fundamental  frequency,  2  the  first 
harmonic,  etc..  Powers  of  ten  in  the  table  are  designated  by  E±xx  which  represents  10*“. 

In  Table  2  changing  the  threshold  over  the  range  l.QxlO^  to  1.0x10'^  has  virtually 
no  impact  on  the  magnitude  of  the  fundamental  and  the  first  three  harmonics,  at  least  to 
four  significant  figures.  The  additional  harmonic  terms  generated  when  the  threshed  is  set 
to  TOxlO**  or  less  appear  to  have  no  impact  on  the  magnitude  of  the  fundamental  and  the 
first  three  hannom'cs.  However,  at  this  point  the  pulse  has  only  propagated  40%  of  the 
distance  to  the  geometric  focal  point  and  very  little  harmonic  generation  has  taken  place. 
A  better  perspective  can  be  obtained  by  examining  the  spectra  at  z  ==  0.9  whidi  is  just 
before  the  geometric  focus.  This  comparison  is  shown  in  Table  3  below.  Again,  to  five 
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sigoificaot  figures,  the  magnitude  of  the  fundamental  and  the  first  four  harmonics  appears 
is  unaffected  by  threshold  values  below  I.OxlO^. 


Table  2.  On-axis  spectra  at  z  =  0.4 


THRESHOLD 

n 

10* 

lO"* 

10  • 

lO** 

10'* 

f 

1.0935E+00 

1.0935E>00 

1.0935E^ 

1.0935E+00 

1.0935E+00 

2 

4.7449E-03 

4.7449E-03 

4.7449E-03 

4.7449E-03 

4.7449E-03 

3 

2.8661E-05 

2.8661E-05 

2.8661E-0S 

2.8661E-05 

2.0568E-07 

4 

5 

6 

7 

2.0568E-07 

2.056aE-07 

2.05e7E-07 

1.5917E-09 

2.0S67E-07 

1.6205E-09 

1.3132E-11 

2.0568E-07 

1.6207E-09 

1.3559E-11 

1.1274E-13 

Table  3.  On-axis  spectra  at 

N 

II 

p 

'O 

THRESHOLD 

n 

10* 

10* 

10* 

10  '* 

10’* 

1 

1.7435E+00 

1.7435E+00 

1.7435E+00 

1.7435E+00 

1.7435E+00 

2 

1.9786E-02 

1 .97e8E-02 

1.9786E-02 

1.9786E-02 

1.9786E-02 

3 

3.2627E-04 

3.2627E-04 

3.2e27E-04 

3.2627E-04 

3.2e27E-04 

4 

6.3104E-06 

6.3097E-06 

6.3O93E-06 

6.3093E-06 

6.3093E-06 

5 

1.1583E-07 

1.3336E-07 

1.3339E-07 

1.3339E-07 

6 

2.8353E-09 

2.9858E-09 

2.9860E-09 

7 

8.8332E'11 

8.9557E-11 

8 

1.6588E-12 

9 

1.8266E-14 

A  last  check  on  the  impact  of  the  threshold  magnitude  is  shown  in  Table  4  which  lists  the 
on-axis  amplitude  ^ctra  40%  beyond  the  geometric  focus  (z  =  1.4).  Again,  for  a  threshold 
between  l.Qxlff*  and  1.0x10'*  there  b  no  significant  impact  on  harmonic  magnitude. 
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Table  4.  On-axis  spectra  at  z  =  1.4 


THRESHOLD 


n 

10-* 

10* 

10* 

10  '• 

io« 

1 

1.2942E+00 

1.2942E4O0 

1.2942E«^00 

1 .2942E+00 

1.2942E^ 

2 

2.t359E-02 

2.1359E-02 

2.1359E-02 

2.13S9E-02 

2.1358E-02 

3 

5.3937E-04 

5.3937E-04 

5.3937E-04 

5.3937E-04 

5.3937E-04 

4 

1.6275E-05 

1.e269E-05 

1.6289E-05 

1.6269E-05 

1.6269E-05 

5 

6 

7 

8 

9 

10 

5.4167E-07 

S.4018E-07 

1.9115E-08 

5.7185E-10 

5.4017E-07 

1.9043E-08 

7.0054e-10 

2.7133E-11 

5.401 7E-07 
1.9043E-08 
6.9950E-10 
2.6488E-11 
1.0397E-12 
2.8673E-15 

5.3  Numerical  Solution  Convergence 

To  assure  that  the  numerical  solution  converges  the  transverse  and  longitudinal  step 
size  are  reduced  until  the  change  in  the  numerical  solution  is  insignificant.  Figure  4  shows 
the  on-axis  normalized  peak  pressure  amplitude  as  a  function  of  the  radial  step  size  used 
in  finite  difference  solution  of  the  KZK  equation. 

In  this  calculation  the  phase  of  the  initial  pressure  is  corrected  using  a  quasi-optical 
approximation.  As  dr  deaeases  from  a  normalized  value  of  2%  of  the  transverse  radius  (dr 
=  0.02)  to  1%  (dr  =  0.01)  the  peak  pressure  profile  changes  dramatically.  However,  as  the 
medi  size  decreases  the  solution  essential^  converges,  at  least  to  a  visual  inqwction.  Tlie 
rate  of  convergence  is  shown  in  Table  5  wdiich  contains  the  value  of  dr,  and  the  magnitude 
of  the  on-axis  pressure  fundamentaL  2”**,  and  3"*  harmonics.  As  dr  ttecreases  the  harmonics 
converge  to  their  asymptotic  values.  By  the  time  that  dr  =  0.004  die  locatkm  d  the 
predicted  focal  point  has  changed  and  the  magnitude  of  each  haimmiic  has  converged  to 
a  point  that  the  jitter  is  on  the  order  of  a  few  parts  per  thousand. 
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GAIN 

1. 


dr 

0.0200 

0.0100 

0.0075 

0.0050 

0.0040 

0.0030 


Geometric 
Focal  Point 


0.004 

0.01 

0.02 


F^re  4.  On-axis  nonnal^d  peak  pressure  convergence 


Table  5.  Solution  convergence  at  predicted  focal  point 


z 

7.4247E-0001 

0.2500E-0001 

8.5439E-0001 

8.5000E-0001 

8.5439E-0001 

8.S439E-0001 


Ann) 

1.7489E*0000 

t.7519E«0000 

1.7532E«0000 

1.7528E40000 

1.7529E40000 

1.7532E^OOOO 


i.aaoiEUooE. 

i.nsoiE-eont^ 

1.8421E-0002V 

1.8312E-0002 

1.8354E-0002 

1.8366E-0002 


An(3 

2.8188E-0004 

2.7882E-0084 

2.8030E-0004 

2.7893E-0004 

2.7822E-0004 

2.78fiOE-0004 
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To  obtain  a  better  approximation  for  the  pressure  in  the  focal  region  an  ampKtude 
shading  can  be  appUed  as  discussed  in  Section  4.  The  calculated  on  axis  peak  pressure 
using  both  phase  and  amplitude  shading  is  compared  in  Figure  5  to  the  calculated  peak 
pressure  using  only  phase  shading.  The  difference  for  the  30°  aperture  is  only  a  few 
percent  in  peak  amplitude  at  the  focal  point  and  the  location  of  the  calculated  focal  point 
does  not  change  as  a  result  of  amplitude  shading. 


Geometric 
Pocel  Point 


Figure  5.  AmpUtude  and  phase  shading  comparison 
6.0  KZK  Line  focus  equation  numerical  predictiems 

This  section  presents  the  predicted  pressure  profile  for  the  experimental 
configuration.  The  inputs  for  the  solution  are  summarized  in  section  6.1.  The  predicted 
peak  pressure  profile  on  and  off  the  symmetry  axis  are  shown  in  Section  6.2.  The  change 
in  the  pressure  time  histoiy  due  to  non-linear  effects  is  shown  in  Section  6.3  when  the  peak 
initial  pressure  across  the  aperture  is  changed  by  an  order  of  magnitude. 
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6.1  Initial  conditions 


The  normalized  initial  peak  pressure  profile  at  the  face  of  the  transducer  is  shown 
in  Figure  6.  The  line  focus  transducer  enclosed  aperture  is  approximately  30”. 

1 

0.8 
0.6 
0.4 
0.2 

j/«  0 

-0.2 
-0.4 
-0.6 
-0.8 
-1 

0.00  0.50  1.00  1.50 

Normalixed  Pressure  30*^  Line  Focub  Transducer 

Figure  6.  Initial  pressure  profile  at  transducer  face. 


The  peak  pressure  at  the  face  of  the  transducer  is  approximately  58.05  kPa.  The  transducer 
radius  of  curvature  is  0.04445  m  with  a  transverse  radius  of  0.G1151  m.  The  input  driving 
frequency  is  1.0667  MHz.  The  transducer  is  in  water  with:  density,  9.98  KgAn’;  velocity, 
1492  m/s;  and  B/A,  5.5.  The  calculations  were  performed  with  the  thresholdiBg  for  the 
floating  boundary  and  the  dynamic  harmonic  allocation  set  to  10^.  The  transverse  step  size 
was  set  to  4.0x10’. 

6.2  Peak  Pressure  Profile  -  58.05  kPa  source 

A  three  dimensional  plot  of  the  pressure  profile  is  shown  in  Figure  7.  The  vertical  (y)  axis 
is  the  peak  pressure  normalized  to  the  initial  peak  pressure  at  the  face  of  the  transducer. 
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The  horizontal  (x)  axis  is  the  range  from  the  transducer  symmetry  axis.  The  z  axis  is  the 
normalized  distance  along  the  direction  of  propagation. 


Figure  7.  Peak  pressure  profile  transverse  to  the  propagation  direction. 

Individual  peak  pressure  profile  segments  are  shown  in  Figure  8  for  as  the  normalized 
distance  increased  from  z  =  0.4  to  z  =  1.3.  This  shows  the  grown  in  the  on-axis  pressure 
as  the  geometric  focus  is  approached  and  the  subsequent  decrease  in  intensity  after  the 
geometric  focus. 
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Figure  8.  Transverse  peak  pressure  profile  vs  normalized  distanoe. 


Since  the  peak  pressure  is  predicted  to  occur  at  z  »  0.86,  Figure  9  Ohistrates  the  grown  in 
peak  pressure  through  the  region  z  =  0.8  through  z  =  0.92. 


C«te  C«la 


G«ia  CaUt 


Figure  9.  Transverse  pressure  profile  around  focal  pc^t 


6.3  Predicted  Impact  of  Increased  Pressure 

When  the  initial  pressure  is  mcreased  the  growth  in  higher  hannom'cs  is  accelerated. 
The  change  in  the  pressure  time  histoiy  for  initial  pressures  of  58.05  kPa  and  508.5  kPa  is 
shown  in  Figure  10.  The  time  domain  wave  shape  is  becoming  asymmetric  i^ich  is 


22 


1 


characteristic  of  the  growth  in  higher  harmonics,  e.g.  5  significant  harmonics  at  58.05  kPa 
and  10  significant  harmonics  at  580.05  kPa.  At  the  same  time  there  is  a  slight  shift  in  the 
location  of  the  predicted  focal,  as  shown  in  Figure  11. 


560.5  kPa 


Figure  10.  Distortion  in  time-domain  waveform  due  to  harmonic  growth. 


kPa 


kPa 


Figure  11.  Shift  in  predicted  focal  point  due  to  increased  pressure. 
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Suaaary  Page  II 


This  section  contains  a  condensed  version  of  the  theory  developed  in  the 
previous  section;  it  was  presented  by  Dr.  Neighbors  at  the  1993  International 
Symposium  on  Nonlinear  Acoustics  in  Norway. 

Also  contained  in  this  section  is  a  condensed  version  of  a  paper  presented 
by  Dr.  Ruf  at  the  1993  Ultrasonics  International  Conference  in  Austria.  It 
details  the  interpretation  of  the  pressure  distribution  within  a  focused  sound 
field  by  making  use  of  the  acousto-optic  Schlieren  image  produced  by  this  field 
and  the  use  of  a  CCD  detector  to  map  this  field.  The  Schlieren  image  evaluation 
and  the  conventional  pinducer  field  mapping  techniques  were  compared  and  the 
field  distortion  due  to  the  presence  of  the  pinducer  probe  were  noted. 


FOCUSED  ULTRASOUND  IMAGING 


Thomas  H.  Neighbors,  Herbert  J.  Ruf,  and  Walter  G.  Mayer 
Ultrasonics  Research  Laboratory,  Physics  Department 
Georgetown  University,  Washington,  D.C.  20057-0995  USA 

ABSTRACT 

Today,  focused  sound  field  measurements  typically  involve  needle  car  membrane 
type  transducers  stepped  across  the  sound  field.  This  time  consuming  process 
is  limited  by  probe  linearity,  effective  aperture  size,  and  experimental 
alignment  In  the  past,  Schlieren  imaging  has  provided  an  effective, 
qualitative  technique  for  examining  a  sound  pressure  field.  This  paper  extends 
the  Schlieren  techniqvie  to  provide  quantitative  measurements  of  the  entire 
focused  sound  pressure  field  without  in-situ  disruptions.  Theoretical  sound 
pressure  field  predictions  for  a  line  focused  system  are  made  based  on  a 
numerical  solution  of  the  Khokhlov-Zabolotskaya-Kuznetsov  parabolic 
equation.  Enhancements  in  computational  effectiveness  of  the  numerical 
solution  are  discussed  and  key  features  of  the  experimental  setup  outlined. 
The  theoretical  predictions  are  compared  with  experimental  measurements 
made  with  the  extended  Schlieren  technique. 


1.  Introduction 

Since  Raman-Nath'  theory  was  published  60  years  ago,  acousto-optic 
techniques  have  been  used  for  non-interactively  measuring  sound  pressure 
fields.  The  diverse  applications  have  ranged  from  the  examination  of  harmonic 
growth^  to  sound  pressure  Held  reconstruction  using  tomography  techniques. 
Schlieren  imaging  as  developed  by  Hiedemann^et  al  has  been  used  extensively 
to  image  the  entire  sound  field  through  photographic  images. 

The  propagating  sound  pressure  for  a  line  focused  radiator  is  predicted 
using  the  Khokhlov-Zabolotskaya-Kuznetsov  (KZK)  parabolic  equation. 
Enhancements  in  the  numerical  solution  are  outlined  and  the  experimental 
setup  for  the  enhanced  Schlieren  technique  summarized.  Predicted  and 
measured  on-  and  off-axis  pressure  proHles  are  then  compared. 


2  Lane  Focused  KZK  Equati(m 

The  KZK  parabolic  equation,  Eq.  (1),  is  used  to  describe  the  propagation 
of  a  directive  fmite  amplitude  sound  beam  including  nonlinearity,  diffraction 
and  attenuation^.  P  is  the  sound  pressure;  z  is  the  on-axis  distance;  t'  is  the 


retarded  time  (t  -  a/Co),  with  the  isentropic  speed  of  sound;  is  the 
transverse  Laplacian;  is  the  fluid  density;  and  p  is  the  coefficient  of 
nonlinearity. 

IL  .  f?vV  -  -J-iL  -  m 

*at'  2  ‘  2p.c*at'» 

The  term  corresponds  to  di^raction  effects.  In  Eq.  (1)  a  well  collimated 
beam  is  assumed,  e.g.  ka  »  1  with  k  =  u/Co  and  a  =  the  source  radius. 

For  a  line  source,  which  considers  only  the  yz  plane  and  ignores  yx  plane 
variations,  the  transverse  Laplacian  can  be  replaced  with  Eq.  (2)  is  the 
line  source  KZK  equation.  It  has  been  simplified  by  replacing  the  retarded 
time  by  a  normalized  retarded  time,  i.e.  t  =  cot';  the  pressure  Pby  a  normalized 
pressure  p  =  P/p^,  where  P(,  is  the  initial  on-axis  peak  pressure;  manipulating 
a  few  terms;  and  using  standard  definitions  of  the  Rayleigh  distance,  z^ 
discontinuity  distance,  1^,  and  the  thermal  viscous  loss  coefficient  a, 

-  cid^  =  (2) 

dodx  42o  ar*  dx*  2/^ 

If  p  is  expanded  as  a  Fourier  series  as  indicated  in  Eq.  (3)  where  is  the 
complex  Fourier  coefficient,  Eq.  (2)  becomes  the  series  of  coupled  partial 
differential  equations  for  the  Fourier  amplitude  coefficients,  shown  in  Eq.  (4). 


. 

'"air- 
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4^0  dr* 


-  acf(/a)V.  =  £/»*  F,.4  ♦  2  £  p^  p/.,  (4) 


Since  the  fundamental  frequency  corresponds  to  n  =  1,  harmonic 
attenuation  is  determined  by  an*.  To  account  for  non-square-law  attenuation 
replace  n*  by  n‘  where  g  is  the  attenuation  rate. 

Since  a  general  analytic  solution  of  the  KZK  equation  has  not  been 
found  to  date,  Eq.  (4)  is  usually  solved  by  numerical  techniques.  The  approach 
used  here  has  been  adapted  from  the  numerical  solution  initially  developed  by 
Aanonsen’,  enhanced  by  Hart  and  Hamilton*,  and  revised  by  Neighbors  and 
Bjorno*.  Computational  efficiency  is  enhanced  by  floating  the  transverse 


boundary  condition  and  dynamically  allocating  harmonic  terms. 

A  traditional  problem  in  Hnite  differenoe  calculations  is  the  calculation's 
outer  mesh  size.  The  simple  approach  shown  in  Figure  1  prevents  the 
inclusion  of  an  artificial  boundary  and  improves  computational  efficiency. 
When  the  pressure  fundamental  frequency  magnitude  at  the  third  point  before 
the  edge  of  the  boundary,  A[n*2],  exceeds  a  pre-specified  threshold  value, 
additional  cells  are  added.  This  causes  the  boundary  to  move  as  the  pressure 
increases  near  the  boundary.  This  approach  prevents  artificial  reflections  at 
the  numerical  boundary  and  the  calculation  time  is  reduced  by  up  to  a  factor 
of  two,  depending  on  the  problem. 

The  second  problem  encountered  in  the  numerical  solution  of  nonlinear 
equations  is  the  large  number  of  harmonic  terms  required  due  to  harmonic 


Figure  1.  'Floating"  transverse  boundary  condition. 

growth.  In  general,  computational  efficiency  is  proportional  to  the  number  of 
harmonics  squared.  A  simple  solution  is  available  consistent  with  the  physics 
of  the  process.  At  the  beginning  of  the  computation,  a  minimum  number  ot 
harmonics  are  included.  At  each  longitudinal  step  in  the  calculation,  the  last 
harmonic's  magnitude  is  checked.  If  it  is  above  a  preset  threshold,  another 
harmonic  is  added.  This  improves  the  computational  efRdency  by  several 
factors  without  any  loss  in  numerical  precision  or  stability. 


3.  Enhanced  Schlieren  Imaging 

Figure  2  illustrates  the  experimental  setup*.  Light  from  a  20  mW,  He- 
Ne  Laser  is  incident  on  a  short  focal  length  lens  which  has  its  focal  point  at 


Filfure  2.  Z-configuration  Schlieren  image  system 

the  spherical  mirror's  focal  point  (4  m).  Thus  the  reflection  from  the  mirror 
is  collimated.  When  the  beam  passes  trough  a  water  tank  containing  the  line 
focus  transducer,  it  is  diffracted  by  the  focused  sound  field.  The  light  is 
focused  by  the  second  spherical  mirror  into  a  diffraction  pattern  on  the  Fourier 
plane.  By  blocking  all  but  the  positive  and  negative  first  diffraction  orders  at 
the  Fourier  plane,  the  sound  field  image  contained  is  incident  on  the  Charge 
Coupled  Device  (CCD)  array  at  the  image  plane.  The  magnitude  of  the 
pressure  to  first  order  is  given  by  the  light  intensity,  1,  =  J|  (w)  where  «  = 
2nyL/X  is  the  Raman-Nath  parameter  with  p  the  maximum  variation  in  the 
index  of  refraction  caused  by  the  pressure,  L  the  width  of  the  sound  beam,  and 
X  the  light  wavelength.  The  CCD  detector  used  in  the  experiment  is  a 
Fairchild  Weston  Systems  model  CCD181,  which  has  2592  active  elements  in 
a  linear  array  with  a  dynamic  range  of  75(X):1  relative  to  rms  noise.  The  line 
focused  transducers  were  constnuled  by  Channel  Industries  of  PZT'4  with  a 
thickness  of 0.625  cm  and  a  radius  of  curvature  of  4.445  cm.  During  the  series 
of  experiments  conducted  three  aperture  angles  of  30^,  6(lP,  and  90^  were  used. 

4.  Comparison  Between  Theixry  And  Bxjieriment 

The  peak  pressure  profile  at  the  face  of  the  30^  aperture  transducer  is 
shown  in  Figure  3  normalized  to  the  on-axis  pressure  at  the  transducer  face 
which  w'as  77.31  kPa.  The  transducer  radius  of  curvature  is  4.445  cm  with  a 
transverse  radius  of  1.151  cm.  The  input  driving  frequency  is  1.1  MHz  and  the 
tank  was  filled  with  water.  The  numerical  values  of  the  density,  isentropic 
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Figure  3.  Initial  pressure  profile  at  transducer  face 
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Figure  4.  On- axis  peak  pressure  experimental  vs  predicted  profiles. 


To  predict  the  on-  and  off-axis  pressure  using  the  line  source  KZK 
approximation,  the  amplitude  and  phase  at  the  transducer  face  were  adjusted 
using  a  quasi  optical  approximation.  The  calculations  were  performed  with  i  he 
threshold  for  the  floating  boundary  and  harmonic  allocation  set  to  10^  and  a 
transverse  step  size  of  4.0x10^.  Figure  4  is  a  comparison  of  the  predicted  and 
measured  on-axis  sound  pressure  for  a  60.19  kPa  and  77.31  kPa  sources.  The 
dashed  line  is  the  theoretical  prediction  usrng  the  line  focus  KZK  equation. 
The  solid  line  is  the  experimental  result  measured  using  the  extended 
Schlieren  technique.  The  difference  between  the  predicted  and  measured  focal 
point  is  approximately  4%  and  the  variation  in  peak  pressure  is  between  10% 
(60.19  kPa  source)  and  1%  (77.31  kPa  source).  In  both  cases  the  quantative 
agreement  between  theory  and  experiment  is  relatively  good. 
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Acousto-Optic  Measurement  of  Pressure  Profiles  in  Ultrasonic 

Fields 


Herbert  J.  Ruf 

Physics  Departnent,  Georgetown  University,  Washington,  DC  20057 


A  Charge-Coupled  Device  (CCD)  detector  is  used  in  a  Schlieren 
imaging  system  to  determine  pressure  profiles  in  an  ultrasonic 
field.  The  ultrasonic  field  studied  is  from  a  1.10  Mhz  line  focus 
transducer  in  water,  operated  in  the  Raman-Nath  regime.  The 
Schlieren  images  are  made  using  only  the  first  diffraction  orders. 
Comparison  of  pressure  profiles  made  with  a  conventional  PVDF  needle 
hydrophone  are  used  to  demonstrate  the  usefulness  of  the  Schlieren 
imaging  technique. 


INTRODUCTION 

Schlieren  imaging  has  been  used  for  decades  to  produce  photographs  of  ultra¬ 
sonic  fields  and  to  investigate  ultrasonic  propagation  through  studies  of  the 
diffraction  patterns  produced.  Schlieren  systems  have  been  utilized  by  Haran 

[1]  to  measure  the  total  acoustic  power  in  an  ultrasonic  beam  and  by  Reibuld 

[2]  to  determine  the  pressure  distribution  in  an  ultrasonic  field.  The  work 
reported  here  demonstrates  a  method  by  which  the  axial  pressure  profile  of  a 
line  focus  transducer  cem  be  measured.  Advantages  of  this  method  lie  in  the 
spatial  resolution  of  the  optical  detector  and  in  the  fact  that  no  mechanical 
scanning  is  necessary. 


THEC«ETICAL  OWISIDERATIOIS 

In  a  Schlieren  imaging  system  a  collimated  beam  of  light  passes  through  a 
region  of  varying  refractive  index.  The  portion  of  the  light  which  is  scat¬ 
tered  by  the  variation  in  refractive  index  forms  an  image  of  the  index  varia¬ 
tions.  In  this  paper  the  index  variations  are  caused  by  the  pressure  gradient 
induced  in  water  by  an  ultrasonic  field.  The  interaction  of  a  progressive, 
sinusoidal  ultrasonic  wave  with  a  light  beam  of  normal  incidence  was  first 
discussed  by  Raman  and  Nath  [3].  The  pressure  gradient  behaves  as  a  phase 
grating  which  produces  a  diffraction  pattern  where  the  intensity  of  light  in 
the  nth  order  is  given  by 

Jj  is  the  nth-order  Bessel  function  and  u  is  the  Raman-Nath  parameter,  which 
is  proportional  to  pressure.  If  only  the  light  from  the  ±nth  diffraction 
orders  is  allowed  to  pass  a  spatial  filter  to  form  an  image,  then  the  light 


intensity  at  each  point  of  the  inage  is  proportional  to  ^  in  eq.  (1). 
Thus,  in  principle,  given  an  image  of  an  ultrasonic  field,  one  can  determine 
the  pressure  field  from  which  this  image  arose. 

Since  v  is  proportional  to  the  pressure  produced  by  an  ultrasonic  trans¬ 
ducer,  and  the  pressure  is  proportional  to  the  excitation  voltage  ai^lied  to 
the  transducer,  V,  the  image  intensity  is  a  function  of  V.  By  varying  V,  and 
recording  the  intensity  of  any  point  in  the  Schlieren  image  as  a  function  of 
V,  data  will  be  acquired  which  will  allow  a  determination  of  the  n  value  at 
this  point,  and  therefore  the  pressure. 


EXPERIMENT 

A  Schlieren  imaging  system  in  a  Z-conf iguration  [4]  consisting  of  a  He-Ne 
laser,  a  spherical  mirror  of  4m  focal  length  as  the  focusing  element,  and  a 
spatial  filter  consisting  of  two  rectangular  slits  which  allows  the  11  dif¬ 
fraction  orders  to  pass,  were  used  for  this  experiment.  The  optical  detector 
used  in  the  image  plane  was  a  2592-element,  Fairchild  linear  CCD  array,  with 
a  lOpm  distance  between  pixel  centers.  The  active  area  of  the  CCD  array  is 
2. 54  cm  and  the  Schlieren  system  was  constructed  so  as  to  allow  a  8.2  cm 
region  to  be  imaged  by  the  CCD,  resulting  in  a  distance  between  pixels  of 
63um- 

The  ultrasonic  sources  were  constructed  with  PZT-4A  line  focus  crystals, 
with  a  radius  of  curvature  of  4.445  cm  and  a  nominal  resonant  frequency  of  330 
kHz,  the  transducers  were  operated  at  the  third  harmonic  at  a  frequency  of 
1.08  MHz.  Two  transducers  were  used,  the  first  with  a  30*  arc  and  the  second 
with  60*  arc.  The  Raman-Nath  theory  assumes  a  uniform  acoustic  field  in  the 
direction  of  light  propagation.  This  is  not  possible  with  a  circular  aperture 
transducer,  so  line  focus  transducers  were  utilized  giving  a  rectangular 
cross-section  to  the  acoustic  field.  The  focused  beam  allowed  for  accurate 
placement  of  the  CCD  array  and  of  the  needle  hydrophones  used  for  a  conven¬ 
tional  pressure  measurement  comparison.  Data  were  acquired  over  a  range  of  10 
to  250  mV  excitation  voltage,  nominally  amplified  by  50  dB,  in  5  mV  steps,  all 
acquisition  was  done  by  computer  with  a  data  acquisition  time  of  0.16  ms  per 
voltage  step,  all  data  were  acquired  within  20  minutes.  Voltage  was  applied 
for  short  durations  at  each  step  to  limit  heat  generation  by  the  transducers 
emd  to  reduce  ultrasonic  streaming. 


ANALYSIS 

A  typical  Intensity  vs.  V  curve  is  shown  in  Fig.  1.  Note  that  the  intensity 
values  are  the  raw  intensity  values  as  measured  from  the  CCD.  Modeling  these 
data  values  with  the  expression, 

I  =  Cjj/lCjV  +  Cj),  (2) 

one  may  utilize  a  curve  fitting  procedure  to  determine  the  constants  Cj,  Cj, 
and  Cj.  From  these  data  the  argument  of  the  Bessel  function,  o,  can  be  calcu¬ 
lated.  Pressure  can  then  be  calculated  from  »  using  the  relationship 

P  =  (0.65w)/L,  (3) 

where  P  is  pressure,  L  is  the  acoustic  beam  width  that  the  light  traverses. 


The  constant  0.65  in  eq.  (3)  conies  from  an  expression  derived  by  Willard  [S] 
which  depends  on  the  frequencies  of  the  light  and  the  sound  beam,  as  well  as 
the  refractive  index  of  the  medium,  water  in  this  case.  The  beam  width  was 
found  by  scanning  the  transducer  along  the  direction  of  light  propagation  with 
a  needle  hydrophone.  The  product  of  the  ratio  of  the  average  to  peak  pres¬ 
sure,  with  the  transducer  crystal  width,  yields  L.  These  beam  scans  are  shown 
in  Figs.  2  and  3  and  demonstrate  that  there  is  no  more  than  2*  of  beam  spread. 


RESULTS  AND  CONCLUSIONS 

Figures  4  and  5  compare  the  pressure  profiles  generated  by  the  above  procedure 
with  profiles  generated  using  a  needle  hydrophone  in  a  scanning  tfink.  Dis¬ 
crepancies  between  the  optical  and  conventional  profiles  can  be  attributed  to 
two  causes,  1)  there  is  a  certain  amount  of  uncertainty  in  the  conventional 
scans  due  to  possible  uncertainty  in  the  calibration  of  the  hydrophones,  and 
2)  Fig.  6  shows  the  effect  of  the  presence  of  the  hydrophone  on  the  acoustic 
field,  resulting  in  elevated  pressure  values.  There  is  also  a  periodic  struc¬ 
ture  visible  in  the  conventional  data,  in  particular  for  the  30*  transducer, 
which  is  caused  by  standing  waves  between  the  transducer  and  the  hydrophone. 

Pressure  profiles  of  line  focus  transducers  were  obtained  using  a  nonin- 
vasive  optical  technique  and  compared  with  pressure  profiles  obtained  with 
conventional  hydrophone  measurements.  A  pressure  artifact,  on  the  order  of 
5%,  caused  by  the  presence  of  the  conventional  hydrophone  was  demonstrated  for 
measurements  made  in  the  focal  region.  Otherwise  overall  agreement  between 
the  two  methods  was  found  to  be  good. 
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Figure  1  Light  intensity  vs.  transducer  excitation  voltage 


Figure  2  Pressure  profile  in  the  direction  of  light  propagation  for  30**  arc  transducer 
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Figure  4  Pressure  profile  inthe  direction  of  sound  propagation  for  30°  arc  transducer 
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Figure  5  Pressure  profile  in  the  direction  of  sound  propagation  for  60°  arc  transducer 
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Figure  6  Pressure  profiles  with  and  without  the  presence  of  a  needle  hydrophone 


Summary  Page  III 


This  section  deals  with  the  various  advantages  acousto-optic  methods  offer 
for  the  investigation  of  sound  fields,  particularly  sound  fields  which  are 
characterized  by  high  amplitude  and  pulsed  regimes.  The  main  thrust  of  the 
present  grant  was  the  investigation  of  focused  ultrasound;  however,  one  has  to 
realize  that  there  are  great  differences  between  focused  continuous  waves  and 
pulsed  waves,  particularly,  as  far  as  the  mathematical  description  of  the  events 
is  concerned. 

Therefore,  an  investigation  was  undertaken  to  examine  the  possibilities 
of  using  an  optical  analysis  of  finite-amplitude  ultrasonic  pulses,  and  the 
feaslbi 1 ity  of  a  spectral  resolution  of  such  pulses,  using  1 ight  diffract  ion  by 
u 1 1  rasound , 

This  section  contains  reprints  of  two  papers  dealing  with  this  topic;  the 
author  of  the  first  paper,  Dr.  John  Wolf,  was  supported  by  the  Grant  at  the  time. 
He  is  now  a  scientist  at  the  Naval  Research  Laboratory  in  Washington,  DC. 
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Spectral  Filtering  of  Finite  Amplitude  Ultrasonic  Pulses 
by  Plates.  An  Investigation  Using  Light  Diffraction 

J.  W.  Wolf,  T.  H.  Neighbors,  III,  W.  G.  Mayer 

Department  of  Physics,  Georgetown  University,  Washington 


Spectral  Filteriiig  of  Finite  AmpHtiide  Ultrasonic 
Putoes  by  Plates.  An  Investigation  Using  Light 
DiffiracthNi 
Saaunary 

Finite  ain{ditude  ultrasonic  pulses,  tefkcted  and  trans¬ 
mitted  from  a  brass  plate,  are  optically  examined 
using  light  diffraction  produced  by  the  ultrasonic  sig¬ 
nal.  Pulse  Fourier  spectra  are  derived  which,  when 
used  for  input  to  a  light  diffraction  model,  provide 
quantitative  agreement  between  experiment  and  the¬ 
ory.  The  spectra  are  used  in  a  propagation  model 
based  on  Burgers’  equation  which  describes  finite  am¬ 
plitude  acoustic  waves  in  a  nonlinear  medium.  Theo¬ 
ries  of  ultrasonic  beam  reflection  and  transmission 
from  plates  immersed  in  a  liquid  are  applied  to  pulsed 
beams  and  incorporated  into  the  pulse  spectrum, 
propagation,  and  light  diffraction  models.  Examples 
are  presented  which  illustrate  selective  filtering  of 
pulse  fundamental  and  harmonic  frequency  compo¬ 
nents  by  Lamb  and  Rayleigh  mode  excitation  in  a 
plate. 

Spektrale  Filtening  von  Ultraschallimpulsen 
endlicher  Amplitude  durch  Flatten,  eine 
Untersuchung  mit  Hilfe  der  Lichtbeugung 
Zusammenfassung 

Es  werden  Ultraschallimpulse  endlicher  Amplitude,  die 
von  einer  Messingplatte  reflektiert  und  transmittiert  wer¬ 
den  mit  Hilfe  der  vom  Ultraschall  hervorgerufenen  Licht¬ 
beugung  optisch  untersucht,  Hierzu  werden  Fourier- 
spektrcn  von  Impulsen  abgeleitet.  die,  als  EingangsgroBe 
fiir  ein  Lichtbeugungsmodell  benutzt  eine  quantitative 
Ubereinstimmung  zwischen  Experiment  und  Theorie  er- 
geben.  Diese  Spektren  werden  in  einem  Ausbreitungsmo- 


dell  verwendet.  das  auf  der  Burgersschen  Gleichung  zur 
Beschreibung  von  Schallwellen  endlicher  Amplitude  in 
einem  nichtlinearen  Medium  beruht.  Die  Theorien  der 
Reflexion  und  Transmission  von  Ultraschallstrahlen  an 
Flatten  in  einer  Fliissigkeit  werden  auf  gepulste  Strahlen 
angewandt  und  in  die  Modelle  fiir  das  Pulsspektrum.  die 
impulsausbreitung  und  die  Lichtbeugung  eingebaut. 
Durch  Beispiele  wird  die  selektive  Filterung  der  Grund- 
frequenz  und  der  Harmonischen  des  Impulses  durch  An- 
regung  von  Lamb-  und  Rayleighwellen  in  der  Platte  de- 
monstriert. 

Etude,  par  ia  diffraction  de  la  iumiere, 
du  filtrage  spectral  d'impulsions  ultrasonores 
d'amplitude  finie  par  des  plaques  solides 
Sommaire 

On  a  examine  optiquement  la  Iumiere  diffractee  par  le 
signal  ultrasonore  reflechi  ou  transmis  par  une  plaque  de 
iaiton  inson^  par  des  impulsions  acoustiques  d'ampli- 
tudes  finies.  Parallelement  on  a  calcule  les  spectres  de 
Fourier  pour  les  impulsions  et  on  les  a  utilises  comme 
signaux  d'entrre  dans  un  modele  de  diffraction  acous- 
tique  de  la  Iumiere.  obtenant  ainsi  un  accord  quantitatif 
entre  calculs  et  experiences.  Le  modele  propagatif  utilise 
repose  sur  I'equation  de  Burgers  decrivant  la  propaga¬ 
tion  d'ondcs  acoustiques  d'amplitude  finie  dans  un  mi¬ 
lieu  non  lineairc.  Puis  on  a  applique  a  des  faisceaux 
puIsM  les  theories  connues  de  la  reflexion  et  de  la  trans¬ 
mission  d'un  faisceau  d'ultrasons  par  une  plaque  immer- 
gw  dans  un  liquide  et  on  les  a  incorporees  dans  des 
modeles  composites  englobant  spectres  d'impulsions, 
propagations  et  diffraction  de  la  Iumiere.  Enfin  on  pre¬ 
sente  des  exemples  illustrant  le  filtrage  selectif  des  com- 
posantes  spectrales  fondamentales  et  harmoniques  des 
impulsions  par  les  modes  de  Lamb  et  de  Rayleigh  excites 
dans  la  plaque. 


Introduction 

Light  diffraction  has  been  used  effectively  as  a  tool  for 
in-situ  measurements  of  ultrasonic  fields  since  the  phe¬ 
nomenon  was  first  observed  by  Lucas  and  Biquard  [1] 
and  Debye  and  Sears  [2],  and  explained  by  the  Raman 
and  Nath  [3]  theory  in  the  1930’s.  By  the  late  1950's, 


Received  30  May  1989, 
accepted  31  December  1989. 

Dr.  J.  W.  Wolf.  T.  H.  Neighbors,  III,  Prof.  W.  G.  Mayer, 
Ultrasonics  Research  Laboratory,  Department  of  Physics. 
Georgetown  University,  Washington,  D.C.  20057  USA. 


Raman-Nath  diffraction  was  successfully  used  to 
probe  continuous  wave  (CW)  finite  amplitude  ultra¬ 
sonic  beams  [4],  Ten  years  later  the  basic  Raman-Nath 
theory  was  generalized  to  account  for  light  sources 
with  Gaussian  intensity  distributions  [5]  and  the  con¬ 
sideration  of  pulsed  rather  than  CW  ultrasound  [6], 
Motivated  by  experimental  observations  in  the  early 
1980’s  [7],  Raman-Nath  diffraction  theory  was  extend¬ 
ed  to  account  for  arbitrary  pulsed  ultrasonic  wave¬ 
forms  [8]. 

Recently  an  experimental  apparatus  was  developed 
using  a  long  optical  path  with  which  one  can  examine 
a  large  variety  of  ultrasonic  pulses  [9].  The  apparatus 
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is  first  used  to  characterize  transducer  frequency  re¬ 
sponse  to  CW  excitation.  The  empirical  results  are 
used  with  a  computer  algorithm  to  construct  a  model 
of  pulse  Fourier  spectra  as  a  function  of  pulse  dura¬ 
tion,  repetition  and  centre  frequencies.  The  technique 
adequately  models  pulses  at  the  transducer  element  in 
the  liquid.  Mayer  and  Neighbors  [10]  combined  the 
pulse  light  diffraction  theory  with  a  model  of  acoustic 
wave  propagation  in  nonlinear  media  [11],  In  their 
paper,  they  suggested  the  use  of  light  diffraction  for 
estimating  the  magnitude  of  harmonic  distortion  in 
finite  amplitude  ultrasonic  pulses.  The  distortion  ap¬ 
pears  as  asymmetries  in  theoretical  pulse  light  diffrac¬ 
tion  pattern  intensity  distributions.  The  optical  sys¬ 
tem  was  used  to  examine  finite  amplitude  pulse 
harmonic  growth  yielding  results  that  agreed  with  the¬ 
oretical  predictions  [12].  These  preliminary  studies 
[9, 12]  have  shown  the  utility  of  the  diffraction  appara¬ 
tus  for  characterizing  ultrasonic  pulses  and  examining 
changes  in  the  pulse  spectrum  due  to  medium  nonlin¬ 
earities. 

The  present  paper  is  concerned  with  the  optical 
analysis  of  pulsed  ultrasonic  beams  before  and  after 
reflection  and  transmission  from  a  solid  plate  im¬ 
mersed  in  the  liquid  medium.  A  satisfactory  theoreti¬ 
cal  description  of  ultrasonic  beam  reflection  from 
plates  was  provided  by  Pitts,  Fiona  and  Mayer  [13] 
where  they  analyzed  the  plane  wave  reflection  coeffi¬ 
cient  and  revealed  the  conditions  necessary  for  an 
ultrasonic  beam  to  generate  vibrational  modes  in  the 
plate. 

Assume  a  monochromatic  ultrasonic  beam  is  inci¬ 
dent  at  some  angle  upon  a  plate  while  its  amplitude 
is  continuously  monitored  at  a  fixed  location  beyond 
the  reflecting  surface  and  along  the  beam  axis.  If  the 
incidence  angle  is  changed  to  a  mode  critical  angle, 
measurement  would  show  a  decrease  in  the  beam  am¬ 
plitude  because  part  of  the  beam  energy  would  no 
longer  be  geometrically  reflected  from  the  plate.  If  the 
incident  beam  was  polychromatic,  such  as  a  train  of 
ultrasonic  pulses  consisting  of  a  discrete  frequency 
spectrum,  a  measurement  of  the  reflected  beam  would 
show  an  amplitude  decrease  in  those  pulse  frequency 
components  that  individually  satisfy  the  conditions  of 
plate  mode  conversion.  Correspondingly,  a  trans¬ 
mitted  pulse  train  will  contain  those  amplitude  terms 
that  originated  from  the  incident  pulse  spectrum  and 
were  removed  from  the  reflected  spectrum  by  mode 
conversion.  The  plate  acts  as  a  frequency  filter.  From 
this  description,  it  is  obvious  that,  for  a  particular 
plate,  the  frequency  spectra  of  the  reflected  and  trans¬ 
mitted  pulse  may  be  significantly  different  from  the 
incident  pulse  Fourier  spectrum.  If  monitored  with  a 
light  beam,  changes  in  pulse  amplitude  spectra  could 
be  directly  measured  from  changes  in  the  intensity 


distribution  of  the  light  diffraction  pattern  produced 
by  these  pulses. 

The  sections  that  follow;  (1 )  describe  the  experimen¬ 
tal  apparatus;  (2)  summarize  the  relevant  aspects  of 
pulse  modeling  and  light  diffraction  by  pulsed  ultra¬ 
sound;  (3)  look  at  a  model  of  high  amplitude  pulse 
propagation  in  nonlinear  media;  (4)  review  the  theory 
of  pulse  reflection  and  transmission  from  plates;  (5) 
demonstrate  the  use  of  the  pulse  propagation  and 
reflection  models;  and  (6)  provide  observations  of 
spectral  filtering  in  experimental  light  diffraction  pat¬ 
terns.  A  more  detailed  description  of  the  experimental 
apparatus,  procedures  and  theoretical  models  used  for 
this  study  can  be  found  in  the  references  [8-13]. 

Experimental  apparatus 

The  light  diffraction  apparatus  can  resolve  light  dif¬ 
fraction  patterns  produced  by  pulses  with  repetition 
frequencies  as  low  as  40  kHz  [9].  The  system  uses 
cylindrical  lenses  which  produce  a  light  beam  of  lin¬ 
ear  cross  section  that  interacts  with  a  slice  of  the  ultra¬ 
sonic  beam  along  its  symmetry  axis.  A  transducer  of 
2.2  cm  diameter  produces  the  ultrasonic  beam  which 
is  intersected  normally  by  the  light  beam  within  a 
water  tank.  The  resulting  diffraction  pattern  is  mea¬ 
sured  by  a  photomultiplier  tube  light  detector.  The 
signal  input  to  the  ultrasonic  transducer  is  driven 
from  the  amplified  output  of  a  pulse  function  genera¬ 
tor.  A  computer  operated  data  acquisition  system 
controls  the  light  detector  and  pulse  function  genera¬ 
tor. 


Pulse  model  and  light  diffraction  theory 

A  model  of  ultrasonic  pulse  spectra  at  the  source  in 
the  liquid  is  based  on  a  measured  ultrasonic  amplitude 
transfer  as  a  function  of  frequency,  ri(/);  for  this 
study  from  a  nominal  1  MHz  transducer  [9].  In  the 
pulse  mode,  the  function  generator  introduces  ideal 
signals  to  the  transducer  which  are  characterized  by 
pulse  centre  and  repetition  frequencies,  /q  and  /p, 
respectively,  and  burst  number,  which  is  related  to 
pulse  duration.  Theoretical  pulse  spectra  are  con¬ 
structed  from  a  discrete  convolution  of  the  electronic 
signal  amplitude  spectrum  with  t,  (/).  The  resulting 
spectrum  is  used  as  input  to  the  light  diffraction 
model  and  is  expressed  as  a  Fourier  sine  expansion 

OD 

v(x,t)  =  v  £  a„sm[rt{a}^t  -  kj,x)  +  (t>„],  (1) 

n  =  0 

where  cUp  =  27t /p,  kp  =  27t/Ap,  Ap  =  the  pulse  wave¬ 
length  in  the  medium  and  v,  the  Raman-Nath  param- 
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eter  [3],  is  proportional  to  the  pulse  peak  amplitude. 
The  a„  and  0„  are,  respectively,  the  amplitude  and 
phase  of  the  Fourier  component  of  the  pulse.  The 
light  intensity  distribution  as  a  function  of  diffraction 
angle  .9  is  given  by  [8] 

X 

f(9)=  S  /„[sin(ksin  9  —  ffifcp)/(/(sin  9  —  mkp)]^, 

m  —  —  X 

(2) 

where: 


L  Z  (Vj 

=  —  X  r„  =  —  X 

...exp{i[ri(m)0,  +  ...  +  r„4)„  +  ...]} 


rj  (m)  =  w  -  2  -  . . .  -  n  r„  -  . . . , 


k  =  2  k//.,  v„  =  va„,  und  J,  is  the  Bessel  function  of 
order  r.  Eq.  (2)  represents  a  diffraction  pattern  which 
exhibits  distinct  maxima  when 


sin9„=  ±///p. 


(3) 


Model  of  nonlinear  pulse  propagation 

The  pulse  propagation  theory  is  based  on  a  modified 
form  of  Burgers’  equation  to  describe  acoustic  wave 
distortion  due  to  medium  nonlinearities.  The  model 
expresses  the  incremental  change  of  the  n"’  harmonic 
component  of  the  pulse  particle  velocity,  f/„,  by  [11] 

U„  (X  +  dx)  =  t4  (X)  +  |(i  ^  cup/c^)  1^  i  j  Uj  U„  _j 

+  (4) 

where:  dx  is  the  incremental  distance,  P  =  \  +(B/2A) 
the  medium  nonlinearity  parameter,  c  the  sound  ve¬ 
locity,  a  the  medium  attenuation  coefficient  and  *  in¬ 
dicates  the  complex  conjugate.  The  use  of  this  model 
will  be  demonstrated  in  a  different  section  of  this  pa¬ 
per. 


Theoretical  pulse  reflection  and  transmission 
from  a  plate 

The  geometry  of  the  pulse-plate  interaction  is  illus¬ 
trated  in  Fig.  1.  Within  a  liquid,  an  ultrasonic  pulse  is 
incident  at  an  angle  0-,  upon  a  solid  plate  of  thickness 
D.  To  describe  the  interaction,  the  monochromatic 
formulism  is  applied  to  pulses  by  multiplying  each  of 
the  incident  pulse  spectral  components  by  the  appro¬ 
priate  reflection  or  transmission  coefficient.  The  coeffi¬ 
cients,  R(k^  and  T(k^,  respectively,  are  complex 
functions  of  k^  =  k^sinGj,  the  wave  vector  compo¬ 


Fig.  I.  Geometry  of  the  pulse-plate  Interaction. 


nent  along  the  liquid-plate  interface,  with  k(,  the  inci¬ 
dent  wave  number  in  the  liquid.  Pitts  et  al.  [13]  define 
the  plane  wave  reflection  coefficient 

/?«£,)  =  JV/[F,FJ,  (5) 

with  the  definitions: 

=  [^s  ~  2  kj]^  [1  —  cos  P]/sin  P 

+  4klx^y.a[l  -  cos Q]/sin Q  -h  i g k* xjx , 
f,  =  [kf  -  2  k^]^  [1  -F  cos  P]/sin  P 

-F  4 k J x, Xj [1  -F  cos Q]/sin Q  —  igk* xjx , 
N  =  [k,^  -  2  k^j*  +  16  ki  xf  x^  -  kf  xj/x^ 

+  8  [k^  —  2  k^]^  x,  Xj  [1  -  cos  P  cos  Q]/[sin  P  sin  Q] . 

=  K  -  sin^  0/1/,  Ka  =  ko  [(Co/Cd)^  -  sin^  ©Js 

k  =  kocos0„  P  =  Dx,i.  Q  =  Dx^  and  q  =  g (liquid)/ 
Q  (solid),  the  ratio  of  the  medium  densities.  The  values 
C5  and  Cd  are  the  shear  and  compressional  wave  veloc¬ 
ities  of  sound  in  the  solid,  respectively.  The  plane  wave 
transmission  coefficient  is  given  [13]  by: 

TikJ  =  {2iekf  Xd/x}  {[kf  -  2k^]VsinP 

+  4klx^xJsinQ}/[F^F/\.  (6) 

With  measured  values  for  c,  c^,  Cd,  g  (liquid)  and 
0  (solid),  the  two  coefficients  reduce  to  functions  of 
/  =  c/2Kko,  D  and  0,.  Since  the  frequency,  /,  and 
plate  thickness,  D,  are  coupled  in  the  expressions, 
f?(kj  and  T(kJ  reduce  to  functions  of  fD  and  0,. 
The  subscripts,  a  and  s,  of  the  functions  F  in  eqs.  (5) 
and  (6)  refer  to  the  antisymmetric  and  symmetric  vi¬ 
brational  modes  of  the  plate  as  defined  in  the  refer¬ 
ences  [14].  When  either  F,  =  0  or  F,  =  0,  the  functions 
F(kJ  and  T{kJ  produce  so-called  pole-zero  pairs 
which  define  the  conditions  giving  rise  to  Lamb  mode 
conversion  of  the  incident  beam  energy  [13]. 
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Fig.  2.  Set  of  velocity  dispersion  curves  for  a  brass  plate  in 
water  with  characteristic  values  listed  in  Tables  I  and  II. 


Fig.  3.  Time  history  and  Fourier  amplitude  spectrum  of  a 
pulse  with  fo  =  2.0  MHz,  =  125.0  kHz,  exponential  rise 
and  decay  time  of  0.85  ps  and  Mq  =  3.0  cm/s:  a)  prior  to  prop¬ 
agation;  b)  after  propagating  50  cm  in  water  with  medium 
values  listed  in  Table  1;  c)  reflected  from  a  brass  plate  with 
D  =  1.69  cm  at  a  Lamb  mode  critical  angle  of  0-,  =  30' ;  and 
d)  transmitted  from  the  plate. 


Table  I.  Measured  characteristics  for  water 
o  =  1.0  g  cm  ’ 

I  c=  1487  9  m  s  (T  =  22  C)  [15] 

B  . 4  =  5.2  [16] 

I  =  T./  ^  x'  =  2.5  X  10“  Nps’  m  [17] 


Table  11.  Measured  characteristics  for  brass. 

a  =  8.5  g  em’ 
c,  =2240.0  m  s  [18] 

1.,  =  4781.0  m  s  [18] 


Fig.  2  shows  a  set  of  velocity  dispiersion  curves  of 
plate  modes  as  a  function  of  critical  angles  of  inci¬ 
dence  and  the  product  /  D,  calculated  with  the  mea¬ 
sured  values  listed  in  Table  II  for  a  brass  plate  im¬ 
mersed  in  water.  The  curves  represent  the  solutions  of 
Fj  =  0  or  =  0  for  real  and  are  labeled  A^,.  S,^.  A; . 
etc.,  to  distinguish  the  antisymmetric  and  symmetric 
plate  modes  of  vibration.  These  curves  serve  as  a  guide 
for  selective  filtering  of  frequencies  in  pulse  spectra. 

A  sample  run 

For  demonstration,  a  theoretical  pulse  train  is  con¬ 
structed  with  an  initial  peak  amplitude  Uq  =  3.0  cm/s, 
/o  =  2.0  MHz,  /p=  125.0  kHz,  and  an  exponential 
rise  and  decay  time  of  0.85  ps.  The  pulse  travels  within 
the  propagation  model  a  distance  of  x  =  50  cm  in  wa¬ 
ter  using  medium  parameters  listed  in  Table  I.  The  top 
two  sets  of  graphs  in  Fig.  3  show  the  normalized  time 
history  and  amplitude  spectrum  of  the  initial  (a)  and 
propagated  (b)  pulse  where  one  sees  a  significant 
amount  of  harmonic  growth. 

The  pulse  is  now  assumed  incident  at  0,  =  30 
upon  a  brass  plate  of  thickness  D  =  1.69  mm.  The 
pulse  centre  frequency  and  plate  thickness  provide  an 
j  D  value  of  3.38  MHz  ■  mm  which  places  the  pulse 
fundamental  frequencies  across  the  plate  mode, 
(see  Fig.  2).  Fig.  3  shows  the  resulting  reflected  (c)  and 
transmitted  (d)  pulse  time  histories  and  amplitude 
spectra.  The  incident  pulse  is  highly  distorted  by  re¬ 
flection  since  several  frequency  terms  have  been  fil¬ 
tered  from  its  spectrum  including  harmonics  filtered  by 
other  plates  modes,  S2,  A2,  etc.  Frequency  terms  from 
the  incident  pulse  spectrum  which  were  filtered  from 
the  reflected  wave  appear  in  the  transmitted  pulse 
spectrum. 

The  dispersion  curves  of  Fig.  2  indicate  that  the  Aq 
and  So  curves  merge  resulting  in  a  large  region  of 
frequency  bandpass  for  the  plate.  Fig.  4  illustrates  the 
reflected  (top)  and  transmitted  (bottom)  pulse  result- 
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Fig.  4.  Time  history  and  Fourier  amplitude  spectrum  of  the 
finite  amplitude  pulse  of  Fig.  3  b  reflected  (top)  and  transmit¬ 
ted  (bottom)  from  a  brass  plate  at  a  Rayleigh  mode  critical 


angle  of  0,  =  45  . 

Source  0 

\ 

Plate 

V 

Source  b 

I'nht  Beam 

Ultrasonic  Bean, 

Fig.  5.  Ultrasonic  beam  configuration  for  the  two  cases  of 
beam  reflection  (Source  a)  and  transmission  (Source  b)  from 
a  plate. 


ing  from  the  theoretical  pulse  of  Fig.  3  incident  upon 
the  brass  plate  at  0;  =  45"  which  excites  a  Rayleigh 
type  plate  mode  [13].  The  reflected  pulse  spectrum  has 
a  greater  percentage  content  of  higher  harmonics  than 
the  incident  pulse.  The  transmitted  pulse  spectrum 
indicates  that  almost  all  higher  harmonics  within  the 
incident  spectrum  have  been  filtered  by  the  plate.  The 
transmitted  pulse  time  history  shows  little  indication, 
except  for  dissipative  losses,  of  finite  amplitude  distor¬ 
tion  and,  essentially,  only  the  initial  pulse  frequencies 
are  present. 

Observation  of  spectral  filtering  by  light  diffraction 

The  actual  orientation  of  the  ultrasonic  source,  plate, 
and  light  beam  in  the  water  tank  for  analysis  of  pulsed 


beams  is  illustrated  in  Fig.  5  for  the  case  of  pulses  (a) 
reflected  and  (b)  transmitted  from  the  plate.  Once  the 
pulse-plate  configuration  is  chosen,  the  data  acquisi¬ 
tion  system  digitizes  the  light  diffraction  pattern  pro¬ 
duced  by  the  acoustooptic  interaction. 

Fig.  6  presents  experimental  light  diffraction  pat¬ 
terns  (left)  and  those  predicted  by  the  combined  pulse 
spectrum,  propagation,  reflection  and  light  diffraction 
models  (right).  The  initial  pulse  was  produced  by  a 
nominal  1  MHz  source  and  was  modeled  using  r,  I  / ) 
and  the  values  —  132.2  kHz,  Jo  =  0,925  MHz,  burst 
number  2,  output  amplitude  Uo  =  3-56cm  s  and 
v  =  1.45.  The  first  light  pattern  set  was  produced  by 
the  finite  amplitude  pulse,  travelling  40  cm  without 
plate  interaction.  The  difference  in  intensity  between 
the  second  set  of  positive  and  negative  satellite  diffrac¬ 
tion  orders,  f+u,  f+i*  and  /+15,  (/=  1.719,  1.851 
and  1.983  MHz),  is  due  to  propagation-induced  har¬ 
monic  generation.  It  is  this  light  intensity  asymmetry 
from  which  one  can  estimate  the  level  of  harmonic 
distortion  in  finite  amplitude  pulses  [10,  12].  For 
the  other  three  light  pattern  sets,  a  brass  plate  was 
mounted  11.5  cm  from  the  source  so  that  ultrasonic 
pulses  are  reflected  from  the  plate  into  the  light  inter¬ 
action  region  while  a  total  travel  distance  of  40  cm 
is  maintained,  (see  Fig.  5  a).  The  plate  thickness, 
£>  =  3.16  mm,  and  the  pulse  centre  frequency, 
Jo  =  925.4  kHz  yield  an  j  D  =  2.923  MHz  •  mm.  The 
light  diffraction  patterns  result  from  incident  angles  of 
©i  =  25',  28  ,  and  32"  where  the  dominant  pulse  spec¬ 
tral  terms  pass  through  the  brass  by  excitation  of  the 
plate’s  A,  Lamb  mode  (see  Fig.  2).  Each  pattern  indi¬ 
cates  that  a  different  spectral  term  has  been  filtered 
from  the  incident  pulse  spiectrum. 

Compare  the  top  set  of  light  patterns  in  Fig.  6  to  the 
bottom  three.  For  0j  =  25",  reflection  introduces  an 
intensity  decrease  in  the  diffraction  orders  which 
corresponds  to  the  n  =  6  spectral  term  of  the  pulse 
(/  =  793.2  kHz).  Looking  toward  diffraction  orders 
corresponding  to  pulse  harmonics,  one  sees  the 
has  decreased  to  below  the  measurement  threshold  of 
the  apparatus  because  a  large  contribution  to  its  in¬ 
tensity  is  l+f,.  The  remaining  orders,  £^,4  and  £+,5 
continue  to  show  asymmetry  with  respect  to  Iq  indi¬ 
cating  harmonics  are  still  present  in  the  pulse  spec¬ 
trum.  The  orders  centred  about  Iq,  that  is  £+1,  de¬ 
crease  since  the  largest  contribution  to  their  intensity 
comes  from  the  pulse  spectral  terms  that  were  re¬ 
moved  by  the  plate. 

Since  it  has  a  positive  slope,  the  A,  mode  passes 
upward  through  the  pulse  spectrum  as  the  incident 
angle  is  increased.  In  the  pattern  set  for  0;  =  28°,  the 
light  order  corresponding  to  the  pulse  centre  frequen¬ 
cy,  £+,(/  =  925.4  kHz),  has  dropped  almost  an  order 
of  magnitude  in  intensity  causing  the  £^1 3.  14  and 
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/  j.  1 5  to  decrease  below  the  apparatus  measurement 
threshold.  The  remaining  /  +  j  diffraction  orders  in  the 
measured  pattern,  are  caused  by  slit  diffraction  within 
the  optical  apparatus.  In  the  set  for  0,  =  32  ,  the  n  =  8 
pulse  spectral  term  has  been  filtered  by  the  plate,  re¬ 
moving  the  /  +  H  ( /  =  1 ,058  M  Hz)  and  the  /  + ,  5  orders 
from  the  diffraction  pattern. 

Fig.  7  compares  the  measured  light  patterns  of 
Fig,  6  produced  by  the  reflected  pulses  (left)  to  mea¬ 
sured  patterns  produced  by  the  incident  pulse  trans¬ 
mitted  through  the  same  plate  at  corresponding 
(right),  (see  Fig.  5  b).  In  each  of  the  three  cases,  the 
patterns  indicate  that  spectral  terms  filtered  by  plate 
reflection  appear  in  the  transmitted  pulse  spectrum. 
The  second  light  pattern  set  shows  the  transmitted 
spectrum  describes  an  almost  monochromatic  ultra¬ 
sonic  beam  where  the  plate  has  transformed  a  pulsed 


Fig.  6.  Measured  (left)  and  modelled  (right)  light  diffraction 
patterns  produced  by  a  finite  amplitude  pulse  with  total 
propagation  x  =  40  cm.  The  top  pattern  set  is  of  the  pulse 
without  a  plate  (a)  and  the  others  are  reflected  from  a  brass 
plate  at  position  x  =  11. 5  cm,  with  &,  =  25^  (b),  28  (c)  and 
32^  (d). 


wave  (interacting  with  the  plate  only  part  of  the  time) 
into  a  quasi-continuous  wave.  Through  mode  conver¬ 
sion,  the  plate  continues  ringing  after  the  incident 
beam  is  not  longer  e.xciting  it.  thereby  broadening  the 
pulse  train  almost  into  a  CW  beam. 

The  apparatus  was  also  used  to  illustrate  selective 
filtering  of  harmonics  generated  by  finite  amplitude 
ultrasound  propagation  as  was  done  theoretically  in 
the  sample  run.  Fig.  8  presents  two  measured  diffrac¬ 
tion  patterns  produced  by  a  transmitted  pulsed  beam, 
originally  incident  on  the  brass  plate  at  a  critical  Ray¬ 
leigh  mode  angle  of  0^  =  45  .  The  resulting  light  pat¬ 
terns  indicate  that  the  fundamental  spectral  compo¬ 
nents  of  the  pulse  were  transmitted  through  the  plate. 
In  the  top  graph  of  Fig.  8,  the  intensity  asymmetry  of 
the  /jij.  /+,4  and  diffraction  orders  indicates 
the  presence  of  harmonics  which  result  from  the  pulse 
propagating  an  additional  28.2  cm  beyond  the  plate  to 
the  light  beam. 

To  demonstrate  harmonic  filtering,  the  brass  plate 
was  positioned  closer  to  the  light  beam  while  keeping 
the  total  pulse  propagation  distance  the  same,  in  the 
bottom  diffraction  pattern  of  Fig.  8.  one  can  see  the 
diffraction  orders  are  symmetrical  with  respiect  to  /y 
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Fig,  7.  Measured  light  diffraction  patterns  produced  by  the 
pulse-plate  configuration  of  Fig.  6.  The  left  patterns  are  of  the 
reflected  pulse  and  the  right  are  of  the  transmitted  pulse. 
0i  =  25  (a).  28  (bland  32  (c). 
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Fig.  8.  Measured  dilTraction  patterns  produced  by  the  trans¬ 
mitted  pulse  originally  incident  upon  the  brass  plate  at 
0,  =  45°.  The  plate  is  located  at  .x=  1 1.5  cm  (top)  and 
X  =  28.5  cm  (bottom). 


indicating  the  plate  has  filtered  the  harmonic  terms  in 
the  pulse  spectrum.  The  plate  has  effectively  removed 
the  harmonic  distortion  of  the  finite  amplitude  pulse. 


Conclusions 

The  experimental  apparatus  has  been  used  to  examine 
noninvasively  finite  amplitude  ultrasonic  pulses  re¬ 
flected  and  transmitted  from  plates  in  a  nonlinear  liq¬ 
uid.  The  theory  describing  reflection  and  transmission 
from  plates  has  been  combined  with  the  pulse  spec¬ 


trum.  propagation  and  light  diffraction  models.  The 
complete  model  predicts  pulsed  ultrasonic  beam  spec¬ 
tral  filtering  by  plates  immersed  in  liquids,  which  was 
then  verified  experimentally  with  the  light  diffraction 
apparatus.  By  examining  its  dispersion  curve,  it  has 
been  shown  that  a  solid  plate  can  be  used  to  filter 
selected  frequencies  of  a  pulse  spectrum,  such  as  indi¬ 
vidual  components  or  whole  regions  of  the  spectrum. 
The  results  demonstrate  the  application  of  light  dif¬ 
fraction  to  analysis  of  ultrasonic  pulses  and  their 
interaction  with  both  the  propagation  medium  and 
solid  materials  within  it.  The  technique  has  the  advan¬ 
tage  of  performing  an  effective  spectral  analysis  of  an 
ultrasonic  beam  with  a  minimum  of  signal  processing 
and  virtually  no  measurable  distortion  of  the  ultra¬ 
sonic  field. 
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Finite-amplitude  pulses  are  examined  acousto-optically  using  a  newly  developed  light- 
diffraction  apparatus.  Based  on  an  optical  analysis  of  ultrasonic  transducer  response  to 
continuous-wave  excitation  at  and  near  the  fundamental  frequency,  pulse  Fourier  spectra  are 
derived  for  input  to  a  light-diffraction  model,  providing  quantitative  agreement  between 
experiment  and  theory.  The  diffraction  theory  predicts  that  a  light-diffraction  pattern 
produced  by  a  harmonically  distorted  acoustic  pulse  train  will  exhibit  asymmetry  in  the 
intensity  distribution  with  respiect  to  the  zero  order.  To  simulate  harmonic  distortion,  pulse 
frequency  spectra  are  used  for  input  to  a  computational  model  that  is  based  on  the  Burgers' 
equation  for  propagation  of  finite-amplitude  acoustic  waves  in  a  nonlinear  medium.  The 
spectrum,  propagation,  and  light-diffraction  models  give  a  complete  description  of  light 
diffraction  by  finite-amplitude  pulses  and  provide  good  agreement  with  experimentally 
obtained  diffraction  patterns. 

PACS  numbers.  43.25.Zx 


INTRODUCTION 

Since  initial  experimental  observations  of  light  diffrac¬ 
tion  by  ultrasonic  waves' '  and  the  subsequent  tlicoictical 
explanation  of  these  observations  by  Raman  and  Nath'  in 
the  193()s,  light  diffraction  has  been  an  effective  tool  for  in 
.situ  measurements  of  ultrasonic  fields.  In  the  late  1950.S,  Ra- 
man-Nath  diffraction  was  successfully  used  for  the  optical 
probing  of  finite-amplitude  continuous-wave  propagation, 
providing  direct  measurements  of  harmonic  components.^ 

In  the  late  1960s  the  basic  Raman-Nath  theory  was  gen¬ 
eralized  to  account  for  light  sources  with  Gaussian  intensity 
distributions'  and  the  consideration  of  pulsed  rather  than 
monochromatic  ultrasound.'’  Initial  experimental  observa¬ 
tions  were  made  of  light-diffraction  patterns  using  pulsed 
ultrasound,^  and  based  on  these  observations,  the  Raman- 
Nath  diffraction  theory  has  been  geneialized  by  Neighbors 
and  Mayer"  to  account  for  arbitrary  ultrasonic  waveforms. 
In  another  paper,  Mayer  and  Neighbors’  proposed  a  tech¬ 
nique  for  estimating  harmonic  generation  in  finite-ampli¬ 
tude  ultrasonic  pulses  by  examining  their  light-diffraction 
patterns. 

Recently,  an  experimental  apparatus  has  been  devel¬ 
oped  using  a  long  optical  path  which  allows  one  to  examine  a 
broad  range  of  ultrasonic  pulse  waveforms.  The  results  ob¬ 
tained  have  validated  the  generalization  of  the  Raman-Nath 
theory  to  arbitrary  pulse  shapes.'”  The  apparatus  first  ana¬ 
lyzes  transducer  response  to  continuous-wave-mode  excita¬ 
tion  (cw)  that  is  used  to  predict  the  structure  of  pulse  Four¬ 
ier  spectra  as  a  function  of  their  duration,  repetition,  and 
center  frequencies.  The  technique  provides  a  quantitative 
agreement  of  the  pulse  diffraction  theory  with  actual  diffrac¬ 
tion  patterns. 

This  paper  presents  a  comparison  of  experimental  dif¬ 
fraction  patterns  produced  by  finite-amplitude  ultrasonic 
pulses  with  the  predictions  of  pulsed  ultrasound  diffraction 
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theory.  The  experimental  observations  provide  an  indirect 
noninteractive  technique  for  observing  growth  and  decay  of 
harmonic  components.  Previous  experimental  observations 
by  Zankel  and  Hiedemann,^  in  examining  the  diffrac  ion 
patterns  produced  by  initially  monochromatic,  finite  ampli¬ 
tude  ultrasound  (for  example,  cw  beams),  showed  that  the 
generation  of  harmonics  in  the  ultrasonic  field  produced  an 
asymmetry  about  the  central  diffraction  order.  The  current 
experimental  observations  show  the  same  trends  for  pulses 
and  demonstrate  the  practical  application  of  light  diffraction 
to  analysis  of  ultrasonic  pulsed  fields  in  nonlinear  media. 

The  sections  that  follow  ( 1 )  outline  the  experimental 
setup;  (2)  describe  the  procedure  used  to  characterize  the 
pulse  prior  to  any  significant  effects  of  propagation  based  on 
transducer  output  characteristics  and  summarize  the  rel¬ 
evant  aspects  of  pulsed  ultrasound  diffraction  theory;  and 
( 3 )  introduce  the  finite-amplitude-pulse  characteristics  into 
this  theory  and  illustrate  the  experimental  observation  of 
harmonic  growth  in  finite-amplitude  pulse  propagation.  A 
detailed  description  of  the  experimental  apparatus  and  the 
pulsed  ultrasound  diffraction  theory  can  be  found  in  the  ref¬ 
erences." 

I.  EXPERIMENTAL  APPARATUS 

Figure  1  is  a  drawing  of  the  light-diffraction  apparatus 
that  uses  He-N’e  laser  light  with  A.  =  632.8  nm  and  accom¬ 
modates  a  light  beam  with  an  object-to-image  plane  distance 
of  7  m  and  width  of  9  cm.  A  cylindrical  lens  diverges  the 
light  beam  in  one  dimension  only,  allowing  all  the  laser  light 
to  pass  through  the  acoustic  beam,  ensuring  proper  normali¬ 
zation  of  a  diffraction  pattern.  The  acousto-optic  interaction 
takes  place  within  a  channel  of  4.5-cm  width,  bounded  by 
glass  windows,  allowing  the  light  beam  to  intersect  normal 
to  the  ultrasonic  beam  within  the  interaction  vessel.  The 
transducer  is  an  air-backed  mounted  PZT  crystal  of  2.5  cm 
in  diameter.  It  is  driven  from  the  output  of  an  rf  power  am¬ 
plifier  that  receives  an  input  signal  from  a  computer-con- 
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FIG.  1.  Schematic  of  experimental  system,  approximately  to  scale. 


trolled  pulse  function  generator.  Acoustic  absorbing  rubber 
is  placed  at  the  end  of  the  channel  to  prevent  reflections, 
permitting  only  progressive  ultrasonic  waves  to  interact 
with  the  light.  Beyond  the  interaction  vessel  a  second  cylin¬ 
drical  lens  collimates  the  diffracted  light  emerging  from  the 
image  plane  into  separate  beams.  The  resolution  of  the  sys¬ 
tem  is  such  that  a  cw  beam  of  1  MHz  produces  a  diffraction 


FIG.  2.  Optically  measured  response  function  of  a  nominal  1-MHz  trans¬ 
ducer,  normalized  to  the  peak  output  amplitude. 


pattern  at  the  screen  with  an  order  separation  of  12  cm. 

Light-intensity  measurements  are  performed  with  a 
photomultiplier  tube  detector  positioned  at  the  screen.  A 
w  omputer-controlled  data-acquisition  system  rotates  a  flat 
mirror,  passing  a  diffraction  pattern  across  a  slit  that  opens 
onto  the  detector  tube  element.  Possible  random  system 
fluctuations  can  influence  the  acquisition  of  data  during  a 
single  pass  of  a  light  pattern.  Therefore,  the  computer  scans  a 
single  pattern  three  to  five  times,  averages,  normalizes  to  the 
peak  light  intensity,  and  stores  the  data  in  a  disk  file.  The  one 
computer  fully  manages  experimental  input,  procedure,  and 
output,  which  consists  of  digitized  diffraction  patterns. 


II.  INITIAL  PULSE  ANALYSIS 

The  optical  system  is  first  used  to  model  ultrasonic 
pulses  prior  to  any  significant  harmonic  distortion  by  non¬ 
linear  propagation. The  construction  of  pulse  spectra  is 
based  on  an  empirical  model  of  transducer  response  to  cw 
excitation.  The  data-acquisition  system  is  programmed  to 
send  monochromatic  signals  to  the  transducer,  which  results 
in  Raman-Nath  diffraction  patterns.  The  experimental  sys¬ 
tem  digitizes  light  patterns  of  m  diffraction  orders  and  con¬ 
verts  their  relative  intensities  /,„  into  values  of  the  Raman- 
Nath  parameter  v  by  the  expression^ 

4  =  [•/,..  (v)]=  (1) 

where  J„,  is  the  Bessel  function  of  order  m.  The  value  of  v  is 
proportional  to  the  acoustic  output  amplitude  of  the  trans¬ 
ducer,  that  is  its  response.  While  keeping  amplitude  of  the 
electronic  signal  input  to  the  transducer  constant,  the  data- 
acquisition  and  control  system  increments  the  frequency  / 
and  determines  v{f)  across  the  nominal  bandwidth  of  the 
transducer.  The  resulting  transducer  amplitude  response  as 
a  normalized  function  of /is  shown  in  Fig.  2. 

The  function  generator  introduces  ideal  electronic  sig¬ 
nals  to  the  transducer  system  that  are  characterized  by  burst 
number  (a  positive  integer  describing  pulse  width),  and 
pulse  center  and  repetition  frequencies /,  and/,,  respective¬ 
ly.  Figure  3(a)  illustrates  the  ideal  pulse  time  history  along 
with  its  Fourier  amplitude  spectrum  for/,  =  0.988  MHz,/, 
=  66.0  kHz,  and  burst  number  5. 

A  model  of  the  pulse  train  at  the  source  within  the  liquid 
is.  constructed  by  a  discrete  convolution  of  the  ideal  pulse 
amplitude  spectrum  with  the  transfer  function  of  Fig.  2. 
Since  light  diffraction  by  cw  ultrasonic  beams  yields  no  in¬ 
formation  relating  to  pulse-phase  spectra,  the  convolved  am¬ 
plitude  spectrum  is  combined  with  the  original  ideal  phase 
spectrum,  which  assumes  that  no  phase  distortion  is  intro¬ 
duced  to  the  amplitude  spectrum  by  the  transducer  system. 
The  resulting  pulse  time  history  and  convolved  amplitude 
spectrum  are  shown  in  Fig.  3(b). 

The  validity  of  the  pulse  model  is  demonstrated  by  using 
it  for  input  to  the  light-diffraction  modeF  and  comparing  the 
predicted  diffraction  pattern  to  the  measured  diffraction 
pattern  produced  by  the  actual  pulse  train  in  the  liquid.  For 
the  theoretical  prediction  by  the  model,  v  is  now  considered 
proportional  to  the  pulse  peak  amplitude.  The  input  for  the 
model  is  expressed  as  a  Fourier  sine  expansion: 
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FIG.  3.  Ideal  pulse  generator  output  time  history  with  amplitude  spectrum 
for =  66.0  kHz,/,  =  0,988  MHz.  burst  5  (top),  and  the  time  history  and 
amplitude  spectrum  after  convolution  with  the  response  function  of  Fig.  2 
(bottom). 


-.10  -CO  I)  ^<1  lU 

Diffraction  Order 


V(JC,/)=V  £  a„  sin [n(co^t  -  k^x)  +  (2) 

n  ^0 

where  =  2n/Xp,  and  =  the  pulse  length  in 

the  medium.  The  a„  and  (t>„  are,  respectively,  the  amplitude 
and  phase  of  the  nth  Fourier  component  of  the  pulse.  The 
farfield  light  intensity  distribution  as  a  function  of  diffrac¬ 
tion  angle  6  is  given  by 


=  i  I,r 

m  =■  ~  x 

with  k  =  Ztt/A, 


' sin(/c  sin  0  —  mk^)  \- 
.  ( /c  sin  0  —  mk  )  )  ' 


=  X  ■■■  X 


I  r,  «  —  «  r,  =  -  oc 


■•■exp{/[/-,(w)<!>,  -t-  •••  +r„<l>„  +  •••]}|  ,  (4) 

v„  =  va„,  and  r,  (m)  =  m  —  Zrj  -  ' ' '  -  —  ■  •  ■  . 

Equation  (3)  represents  a  set  of  discrete  diffraction  or¬ 
ders  each  with  an  individual  distribution  described  by  the 
[sin  X /X]  ^  term  which  is  maximum  when 

sin  =  ±mA  /X^.  (5) 

For  small  pulse  amplitudes,  Eq.  (4)  reduces  to  /g  ~  1,  and 

(6) 

Equations  (5)  and  (6)  relate  the  dominant  terms  of  the 
pulse  amplitude  spectrum  directly  to  individual  diffraction 
orders.  For  values  of  v  <  2,  Eq.  ( 6 )  provides  a  direct  estimate 
of  a„  from  the  measured  intensity  I„ .  Figure  4  shows  the 
digitized  pattern  and  resulting  prediction  of  the  pulse  spec¬ 
trum  and  light-diffraction  models  where  good  agreement  is 
seen  in  all  parts  of  the  light  pattern. 


FIG-  4.  Measured  (top)  and  modeled  (bottom)  Ught-dilfraction  patterns 
produced  by  the  pulsed  beam  of  Fig.  3  with  v  =  1.60. 


III.  NONLINEAR  PULSE  PROPAGATION 

Mayer  and  Neighbors’*  calculate  nonlinear  propagation 
of  pulses  with  a  model  presented  by  Haran  and  Cook." 
Based  on  Burgers’  equation,  the  model  describes  wave  pro¬ 
file  distortion  as  the  pulse  propagates  along  the  x  direction, 
including  arbitrary  frequency-dependent  attenuation  and 
harmonic  generation.  The  one-dimensional  Burgers’  equa¬ 
tion  expresses  the  differential  change  of  the  nth  harmonic 
component  of  the  particle  velocity  U: 


where  c  is  the  small  signal-sound  velocity,  t  =  t  —  (x/c)  is 
the  retarded  time,  a  the  medium  attenuation  coefficient,  and 
{B /2A)  is  the  nonlinearity  parameter.  Expanding 
(/  in  a  Taylor  series  and  combining  terms  up  to  first  order 
with  Eq.  (7),  Cook  and  Haran  derive 


U„  (x  -I-  dx)  =  U„  (x)  -f- 


+  i  nUjUf_„\-an^U„\dx.  (8) 

j=n+j  / 

Equation  (8)  gives  the  incremental  change  of  an  acoustic 
wave’s  spectral  components,  where  dx  is  the  incremental 
distance  and  ♦  refers  to  the  complex  conjugate. 

As  input,  the  propagation  model  is  given  a  pulse  spec¬ 
trum,  modeled  with  an  exponential  rise  and  decay,  along 
with  a  set  of  medium  constants  and  initial  conditions.  The 
pulse  has  fg  =  2.0  MHz,/^  =  125.0  kHz,  and  a  rise  and  de- 
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TABLE  II.  Measured  characteristics  for  water 


FIG.  5.  A  finite  amplitude  pulse  of/i  =2MHz,y^  =  125  kHz,  and  an  expo¬ 
nential  rise  and  decay  time  of  0.85  /rs.  Medium  parameters  are  a  =  a'f', 
a'  =  1 .9  X  10  '''  Np  sVm,  B  //I  =  5.5,  c  =  1492  m/s,  and  initial  amplitude 
Uu  =)  cm/s. 


cay  time  of  0.85  /zs.  The  medium  constants  used  for  this 
example  are  5  7^4  =  5.5,  c  =  1492  m/s,  and  a  quadratic  fre¬ 
quency  dependent  attenuation,  a  =  a'/^  where  a' 
—  1 .9  X  10  “  Np  s^/m.  The  initial  condition  consists  of  a  3 
cm/s  peak  particle  velocity  Uq  .  Figure  5  illustrates  the  model 
output,  with  the  pulse  time  history  and  amplitude  spectrum, 
normalized  to  Uq,  at  four  points  of  propagation.  This  illus¬ 
trates  the  expected  harmonic  generation  and  the  associated 
shock-front  formation  in  each  cycle  of  the  pulse.  Also  indi¬ 
cated  is  the  overall  decrease  of  the  pulse  amplitude  by  at¬ 
tenuation. 

To  compare  the  model  to  an  actual  case  of  pulse  propa¬ 
gation,  an  initial  particle  velocity  must  be  assigned  to  the 


TABLE  I.  Peak  amplitude  values  for  the  pulse  of  Fig.  3  extrapolated  from 
known  parameters  at  x  =  4.5  cm  using  Eq.  (9). 


X  (cm) 

<t>  (mV) 

V 

Uo  (cm/s) 

4.5 

81.6 

1.60 

3.12 

10.0 

69.6 

1.35 

2.25 

20.0 

76.4 

1.35 

2.46 

30.0 

88.4 

1.35 

2.85 

40.0 

101.9 

1.35 

3.29 

c=  1487.9  m/s  {T=2TCy 
a  =  a'f',  a' =  2.5x10  '■*  Np  s’/m' 

■“H.  Medwin,  i  Acousl.  Soc.  Am  58.  1318  ( 1975) 

’’H  A  Kashkooli.  P.  J.  Dolan,  Jr  ,  and  C.  W.  Smith,  J.  Acoust.  Soc.  Am  82. 
2086(1987). 

'J.  M.  M.  Pinkerton.  Nature  160.  128  ( 1947). 


pulse.  Zankel  and  Hiedemann  reference  an  expression'’  that 
relates  m„  to  v  and  the  width  of  the  ultrasonic  beam,  D.  From 
Schlieren  images,  a  D  is  estimated  and  provides  w,,  =  3.12 
cm/s  for  the  pulse  shown  in  Fig.  3.  The  optical  beam,  of9-cm 
width,  is  assumed  to  sample  the  ultrasonic  field  at  x  =  4.5 
cm,  the  distance  from  the  transducer  face  to  the  midpoint  of 
the  sound  path  through  the  light-ultrasound  interaction  re¬ 
gion. 

The  Schlieren  method  of  estimating  D  is  acceptable  near 
the  ultrasonic  source,  but  elsewhere  is  not  accurate  due  to 
beam  divergence.  However,  the  beamwidth  is  considered  by 
each  value  of  v.  Also,  u„  is  linearly  related  to  the  input  vol¬ 
tage  <1>  to  the  transducer  system.  These  two  points  are  sum¬ 
marized  in  the  relation 

w„=r«l>v,  (9) 

where  T  is  an  experimental  constant  determined  from  initial 
conditions. 

As  an  example,  the  pulse  of  Fig.  3  is  propagated  while 
keeping  v  constant  for  x  >  4.5  cm.  Appropriate  values  of  w„ 
are  extrapolated  from  initial  pulse  parameters  using  Eq.  (9) 
and  listed  in  Table  I.  The  medium  constants  used  for  water 
are  given  in  Table  II.  The  pulse  spectrum  and  the  values  of 
Tables  I  and  II  are  used  for  input  to  the  propagation  model. 
The  predicted  pulse  amplitude  spectrum  shown  in  Fig.  6 
illustrates  the  expected  harmonic  growth  at  10,  20,  30,  and 
40  cm  from  the  transducer. 

The  level  of  harmonic  distortion  shown  in  Fig.  6  is  small 
compared  to  the  example  of  Fig.  5  since  the  total  propaga¬ 
tion  distance  is  smaller.  This  low  level  of  distortion  presents 
no  problem  for  this  investigation  due  to  the  diffraction  pat¬ 
tern’s  extreme  sensitivity  to  pulse  shape.  The  propagated 
pulse  spectra  are  used  for  input  to  the  diffraction  model  for 
comparison  to  actual  light  patterns. 

Experimental  pulse  propagation  is  performed  with  the 
optical  system  by  positioning  the  transducer,  within  the  in¬ 
teraction  vessel,  at  selected  distances  away  from  the  axis  of 
the  collimated  light  beam.  Once  the  transducer  is  fixed  at  the 
desired  distance,  the  data-acquisition  system  performs  the 
experiment  and  stores  the  resulting  digitized  diffraction  pat¬ 
tern.  The  measured  patterns  and  the  diffraction  model  out¬ 
put  are  shown  in  Fig.  7. 

Analogous  to  the  results  of  Zankel  and  Hiedemann"'  for 
cw  ultrasonic  beams,  the  asymmetry  of  the  diffraction  orders 
about  7o  can  be  seen  in  both  the  measured  and  modeled  pulse 
light  patterns.  The  growth  of  the  first  set  of  pulse  harmonics 
is  indicated  in  the  diffraction  orders  about  I 2  mhz  (^zs  ). 
which,  due  to  a  phase  shift,  appears  as  a  decrease  in  the 
diffraction  orders  about  /_  j  mhz  (^-  28 ) ■ 
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0.5 


lO  cm 


TABLE  111.  Peak  amplitude  values  for  a  pulse  of  —  122.0  kHz. 
ft)  ~  0.987  MHz,  and  burst  number  3. 


X  (cm) 

4>  (mV) 

y 

Uu  (cm/s) 

40.0 

73.4 

0.81 

1.80 

40.0 

100.0 

1.11 

2.46 

40.0 

126.7 

1.41 

3.12 

As  another  example,  a  pulse  of  different  spectral  content 
is  propagated  and  analyzed  optically.  The  pulse  is  sampled 
repeatedly  at  x  =  40  cm;  however,  in  this  case  the  initial 
amplitude  u„  is  varied.  Since  the  same  ultrasonic  source  is 
used,  a  pulse  spectrum  is  generated  using  the  response  curve 
of  Fig.  2,  but  with  values =  122.0  kHz, =  0.987  MHz, 
and  burst  number  3.  The  constants  of  Table  II  are  used  with 
values  for  and  v  summarized  in  Table  III. 

The  combined  output  of  the  propagation  and  diffraction 
models  are  compared  to  measured  light  patterns  in  Fig  8. 
The  harmonic  distortion  increases  with  initial  pulse  ampli¬ 
tude  as  predicted  by  theory.  This  set  of  diffraction  patterns 
illustrates  the  decrease  of  the  original  amplitude  terms  due 
to  the  transfer  of  their  energy  to  the  first  set  of  pulse  harmon¬ 
ics.  Both  examples  illustrate  good  agreement  between  mea¬ 
surement  and  theory,  as  well  as  the  sensitivity  of  optical 
analysis  to  even  subtle  harmonic  distortions  of  pulsed  ultra¬ 
sonic  fields. 

IV.  CONCLUSIONS 

The  results  presented  here  demonstrate  the  applicability 
of  light  diffraction  to  an  analysis  of  ultrasonic  pulses  in  a 
nonlinear  medium.  With  the  pulse  modeling  provided  by 
preliminary  diffraction  studies,  a  pulse  spectrum  has  been 
constructed  for  input  to  the  propagation  model  of  Haran  and 
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Frequency  (MHz)  Diffraction  Order 


FIG.  7.  Measured  (left)  and  modeled  (right)  light-diffraction  patterns  pro-  FIG.  8.  Measured  (left)  and  modeled  (right)  light-diffraction  patterns  for 

duced  by  the  pulsed  beam  spectra  of  Fig.  6  each  with  v  =  1,35.  medium  constants  and  pulse  values  as  listed  in  Tables  II  and  III. 
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Cook."  The  predictions  of  the  propagation  model,  when 
used  with  the  diffraction  algorithm,  show  good  agreement 
with  light  patterns  measured  by  the  experimental  apparatus. 
Although  very  subtle  for  small  propagation  distances,  pulse 
harmonic  distortion  produces  an  observable  asymmetry  in 
light-diffraction  patterns.  These  results  support  the  sugges¬ 
tion  of  Mayer  and  Neighbors'*  that  an  estimate  of  pulse  har¬ 
monic  growth  can  be  attained  by  a  study  of  light-diffraction 
patterns. 
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